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ABSTRACT 

We use a sample of ~16,000 non-emission line galaxies from the SDSS to investigate the physical 
parameters underlying the well-known color-magnitude and color-cr relations. Galaxies are sorted in 
terms of velocity dispersions (cr), luminosity (L), and color, and their spectra are stacked to obtain 
very high S/N mean spectra for stellar population analysis. This allows us to map mean luminosity- 
weighted ages, [Fe/H], [Mg/H], and [Mg/Fe] in cr-L-color space. Our first result is that there are 
many different red sequences, with age, [Fe/H], [Mg/H], and [Mg/Fe] showing different amounts of 
slope and scatter when plotted versus cr, L, or color. These behaviors are explained if the star 
formation histories of the galaxies populate a two-dimensional parameter space. One parameter is the 
previously well-known increase in age, [Fe/H], [Mg/H], and [Mg/Fe] with a. In addition to this, we 
find systematic variations at fixed ct, such that more luminous galaxies are younger, more Fe-rich, but 
have lower [Mg/Fe] than their fainter counterparts. The main a trends support a paradigm in which 
more massive galaxies form their stars more rapidly and at earlier times than less massive galaxies. 
The trends at fixed a are consistent with scatter in the duration of star formation for galaxies at a 
given <j. The co-variation of stellar population properties and L residuals at fixed a that we present 
here has a number of implications: it explains the differing behavior of stellar population indicators 
when investigated versus a as compared to i, and it reveals that L is not as efficient as a for indicating 
galaxy "size" in stellar population studies. 

Subject headings: galaxies: abundances, galaxies: elliptical and lenticular 



1. INTRODUCTION 

Early type galaxies are observed to obey many scaling 
relations between their structural and spectral proper- 
ties. Early work identified a number of one-dimensional 
relations, including a color-magnitude relation (Faber 
1973; Sandage & Visvanathan 1978; Bower et al. 1992), 
color-line strength relations (Faber 1973), the Faber- 
Jackson relation between galaxy luminosity (L) and ve- 
locity dispersion (cr) (Faber & Jackson 1976), variations 
in galaxy mass-to- light ratio (M / L) versus L (Tinsley 
1981; Faber et al. 1987), and various correlations of cr and 
L with galaxy effective radius {Re) and effective surface 
brightness (/Lte) (Kormendy 1985), or with galaxy core 
radius (rc) and central surface brightness (Iq) (Lauer 
1985). It is clear however that early type galaxies ac- 
tually comprise at least a two-parameter population in 
terms of their structure, as represented by the Funda- 
mental Plane (FP) relation (Djorgovski & Davis 1987; 
Dressier et al. 1987). Projections of the FP appear nar- 
row in some orientations, leading to the seemingly one- 
dimensional relations listed above. 

In comparison to structure, our understanding of the 
stellar population scaling laws is less well developed and 
has primarily focused on one-dimensional relations, de- 
spite the evidence that early type galaxies form a two- 
parameter population, at least structurally. The stellar 
populations of early type galaxies are known to vary with 
galaxy "size" , where size is typically parameterized as 
either cr or L, but sometimes also as stellar mass (M*). 
Such relations include the color-magnitude and color-cr 
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relations (where it is assumed that the galaxy color is 
indicative of its stellar population) as well as absorption 
line strength-cr relations (e.g., Burstein et al. 1984; Ben- 
der et al. 1993; Trager et al. 1998; CoUess et al. 1999; 
Kuntschner et al. 2001; Bernardi et al. 2003b), which 
give potentially more detailed information about the stel- 
lar populations than do broad-band colors. In general, 
colors are observed to redden and metal line strengths to 
increase with increasing a or L. 

The question of whether these variations are due pri- 
marily to variations in metallicity, or whether age varia- 
tions may also play a role was first posed by Faber (1973) 
and is still debated in the literature. The age-metallicity 
degeneracy is such that differences in population age 
and differences in population metallicity affect colors and 
metal absorption lines in nearly indistinguishable ways. 
This can be loosely quantified by the "3/2" rule, that a 
fractional change in metallicity is roughly equivalent to 
a 1.5 times larger fractional change in age in terms of its 
effect on colors and metal lines (Worthey 1994). Many 
authors (e.g.. Bower et al. 1992; Kodama & Arimoto 
1997) have since claimed that the color-magnitude re- 
lation is driven primarily by metallicity variations rather 
than age variations, although Faber et al. (1995) argued 
that age effects were also important. Recently, Gallazzi 
et al. (2006) showed that galaxy metallicities increase 
by 100% along the red sequence, while ages increase by 
only 50%. Applying the 3/2 rule, metallicity variations 
should therefore account for ~ 75% of the observed color 
variation along the color-magnitude relation, while age 
variation provides the remaining 25%. 

Other recent work studying stellar population trends 
with cr (rather than L) have demonstrated age variations 
with cr that should make larger contributions to the color 



2 



Graves et al. 



variation. Applying the 3/2 rule to these studies, age 
appears to contribute from ~ 33% (Thomas ct al. 2005; 
Graves et al. 2007) to ~ 40% (Trager et al. 2000; Nelan 
et al. 2005; Smith et al. 2007a) of the observed color 
variation along the color-cr relation. 

In this paper, we tackle these subtle differences head- 
on, comparing the color-magnitude and color-cr relations 
directly. We explore in detail the related stellar popu- 
lation variations, focusing on age, metallicity, and abun- 
dance ratio variations with a and separately with L. We 
then look at how stellar populations correlate with resid- 
uals from the mean a-L relation. This analysis explains 
how these systematic residual relations conspire to pro- 
duce the similar yet subtly and confusingly different ob- 
served trends with a and L. 

This is the first in a series of papers presenting a de- 
tailed exploration of stellar population properties of qui- 
escent (non-star forming) galaxies, as well as the depen- 
dence of these properties on the structural parameters 
of the galaxies. The ultimate goal is to match up infor- 
mation about the star formation histories of quiescent 
galaxies (derived from stellar population analysis) to the 
mass assembly and structural evolution of these galaxies 
(as inferred from galaxy structure and morphology). We 
use a sample of '^16,000 quiescent galaxies from the Sloan 
Digital Sky Survey (SDSS), stacking dozens or hundreds 
of similar galaxies to get very high signal-to-noise [S/N) 
spectra for detailed and accurate stellar population work. 
For the purpose of this analysis, we differentiate between 
"global variables" , which we use to identify similar galax- 
ies for stacking, and "stellar population variables" , which 
are the mean ages, metal abundances, and abundance ra- 
tios derived from the stacked spectra. 

This first paper focuses on the color-magnitude and 
color-cr relations of quiescent galaxies; the global vari- 
ables considered are therefore cr, i, and color^. The goal 
is to understand where these relations come from, what 
causes scatter around the relations, and how and why the 
color-magnitude and color-cr relations differ from one an- 
other. The color-magnitude relation is a natural place to 
begin this series of papers, as it is the most accepted 
and most widely used of the early type galaxy scaling 
laws, to the extent that it is even sometimes used to 
identify galaxy clusters in photometric surveys where no 
redshifts are available (e.g., Gladders & Yee 2000). The 
stellar population variables presented here include mean 
stellar age, total Fe abundance ([Fe/H]), total Mg abun- 
dance ([Mg/H]), and the relative abundance ratio of Mg 
and Fe ([Mg/Fe]). 

Much work has already been done on these topics and 
not all of the results presented in this paper are new. 
However, the origin of the color-magnitude relation is 
fundamental to understanding the formation and evolu- 
tion of early type galaxies, yet as discussed above, there 
is considerable confusion in the literature as to how the 
various scaling laws fit together. By directly investigat- 
ing the various relations between ct, L, color, and stel- 
lar population properties simultaneously, we show that 
the various stellar population variables correlate quite 

* It is worth keeping in mind that, although L and color are 
treated as "global" variables in this work, they (unlike a) depend 
significantly on the stellar population properties of the galaxy and 
are thus not truly independent structural parameters. We call them 
"global" (along with a) for linguistic simplicity. 



strongly with some of the global variables but only very 
weakly with other global variables. In particular, a and 
L are not interchangeable. For instance, age correlates 
strongly with a but rather weakly with L, while [Fe/H] 
shows the opposite behavior. 

Briefly looking forward to the subsequent papers in 
this series, Paper II maps stellar population variations 
throughout the Fundamental Plane (FP). In Paper II we 
show that the stellar populations of galaxies on the FP 
vary only with ct; the second dimension of the face-on FP 
appears to be uncorrelated with the galaxy star forma- 
tion histories. However, variations are observed through 
the thickness of the FP. Paper HI demonstrates that the 
2-D family of quiescent galaxy stellar populations can 
be mapped onto a 2-D cross-section through the FP and 
presents a toy model of galaxy star formation histories. 
Finally, Paper IV compares various stellar (M^/L) and 
dynamical {Mdyn/L) mass-to- light ratios for the galaxy 
sample. Paper IV demonstrates that the variations in 
M*/ L due to stellar population effects are inadequate to 
explain the observed variation in Mdyn/L, both along the 
FP and perpendicular to the plane. 

In §2, we describe the sample of galaxies used in this 
analysis. Section 3 reviews known properties of the color 
magnitude relation for early type galaxies and the dis- 
tribution of these galaxies in the three dimensional <j-L- 
color space as a motivation for the subsequent analysis. 
In §4, we discuss the process by which we deconstruct the 
color-magnitude relation and measure stellar population 
properties across the red sequence. Section 5 presents 
the results of the analysis and maps stellar population 
ages, [Fe/H], [Mg/H] and a-enhancement in the three- 
dimensional CT-L-color space. This is followed by a brief 
discussion of the results in §6, and a summary of our 
conclusions in §7. 

2. THE QUIESCENT GALAXY SAMPLE 

The goal of this analysis is to achieve a fundamen- 
tal understanding of the correlations between ct, L, and 
color for typical early type galaxies and their relation to 
stellar population properties. However, from the very 
outset we face challenges in defining a reasonable sample 
with which to explore these relations. Although there 
is a broad correspondance between galaxies with early 
type morphologies, galaxies that lie on the red sequence, 
and galaxies that are not actively forming stars, not all 
early type galaxies are on the red sequence, not all red 
sequence galaxies have early type morphology, and some 
early type galaxies are still actively forming stars. 

Schawinski et al. (2007) identify a sample of morpho- 
logically selected SDSS early type galaxies and classify 
their emission line properties using the emission line ra- 
tio diagrams pioneered by Baldwin et al. (1981, here- 
after BPT). Their early type sample consists of 81.5% 
quiescent galaxies (i.e., no > 3ct emission detections in 
Ha, H/3, [Oiii]A3727, or [Nii]A6583), 1.5% Seyferts, 5.7% 
low ionization nuclear emission- line regions (LINERs), 
and 11.2% star- forming or transition objects (which host 
both star formation and AGN activity). The major- 
ity of their star-forming galaxies, transition objects, and 
Seyferts are not on the red sequence; clearly it would not 
make sense to include such objects in an analysis of the 
color-magnitude and color-CT relations for passive galax- 
ies. Schawinski et al. (2007) require > 3ct detections in 
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all four emission lines listed above and therefore likely 
identify as "quiescent" numerous galaxies with low-level 
emission that may be detected in the strongest emission 
line but is below the 3a threshold in weaker lines such as 
H/3 or [Nil]. 

In a parallel study, Yan et al. (2006) classify a volume- 
limited color-selected sample of SDSS red sequence 
galaxies based on their emission line ratios, also using 
BPT diagrams. They identify as quiescent only galax- 
ies that have no emission lines detected at the 3a level. 
In practice, the strong Ha and [Oii]A3727 emission lines 
drive this definition. They find that only 47.8% of red 
sequence galaxies are quiescent, with a further 28.8% ex- 
hibiting LINER-like emission and 23.4% likely hosting 
Seyfert activity, star formation, or a combination of the 
two. Their red sequence galaxies almost certainly include 
dust-roddencd star- forming galaxies, who dust-corrected 
colors would lie off of the red sequence. 

In light of these studies, we have chosen to study the 
correlations between a, L, color, and stellar populations 
for a sample of quiescent galaxies. Excluding galaxies 
with active star formation removes potential contami- 
nation from dust-reddened star-forming objects whose 
colors do not reflect their underlying stellar populations. 
The quiescent criterion also excludes Seyfcrts and LIN- 
ERs. Many Seyferts have early type morphologies but 
typically have substantially younger stellar populations 
(Kauffmann ct al. 2003), to the extent that they mostly 
do not fall on the red sequence (Schawinski et al. 2007) 
and are not relevant to the typical early type galaxy 
color-magnitude and color-cr relations. Furthermore, a 
strong point source may contaminate global galaxy color 
determinations . 

A substantially stronger case could be made for includ- 
ing galaxies with LINER emission because they are com- 
mon on the red sequence and typically reside in early type 
hosts with red colors. In Graves et al. (2007), we anal- 
ysed the stellar populations of non-star-forming red se- 
quence galaxies, comparing quiescent galaxies with those 
hosting LINERs. We found that the LINER hosts had 
systematically younger ages, but the same metallicities 
and abundance patterns as quiescent red sequence galax- 
ies at the same a. They also have similar or slightly 
redder colors than their riuicscent counterpart, suggest- 
ing that the color effects of their younger populations are 
offset by larger mean dust content. Excluding them from 
the sample in this analysis will remove many of the most 
recent arrivals on the red sequence. 

However, we have; dcicided to limit the analysis in this 
work to quiescent galaxies only, saving the LINERs for 
a future separate analysis. It can be difficult to distin- 
guish between LINERs and transition objects (Schawin- 
ski et al. 2007). There is no real consensus in the lit- 
erature as to the true nature of LINER sources, which 
may well represent a heterogeneous populations. Fur- 
thermore, there is some indication that the LINER sam- 
ple of Graves et al. (2007) in fact contains two sepa- 
rate sub-populations. A preliminary analysis suggests 
that these sub-populations have different stellar popula- 
tion properties from one another and from the quiescent 
galaxy population, as well as different luminosity and 
mass distributions at fixed a. As such, the LINER pop- 
ulation warrant a separate and detailed analysis. 

The sample discussed in this paper is therefore a sam- 



ple of quiescent galaxies, chosen spectroscopically to be 
free of emission lines of any kind. It is not a represen- 
tative sample of "early type galaxies" or "red sequence 
galaxies", but instead provides a window onto the stel- 
lar population properties and star formation histories of 
galaxies that are no longer forming stars. The reader 
should bear in mind that the exclusion of LINERs from 
the sample may remove many galaxies whose star forma- 
tion ended very recently and that the mean ages of the 
galaxies in this analysis are therefore biased toward older 
values. If the stellar population trends of the LINER host 
galaxies follow the same behavior as quiescent galaxies, 
the relative properties presented here should hold for the 
non-star-forming galaxy population as a whole even if 
the zeropoint of the age scale is uncertain. However, the 
properties of the LINER population must be explored in 
future work to verify this possibility. 

2.1. SDSS Data 

The galaxies used in this analysis are drawn from 
the Sloan Digital Sky Survey (SDSS, York et al. 2000), 
which was conducted with a dedicated 2.5m telescope at 
Apache Point Observatory (Gunn et al. 2006). Images 
are taken in drift-scanning mode with a mosaic CCD 
camera (Gunn et al. 1998), processed (Lupton et al. 2001; 
Stoughton et al. 2002; Pier et al. 2003), and calibrated 
(Hogg et al. 2001; Smith ct al. 2002; Ivezic et al. 2004; 
Tucker et al. 2006). All galaxies used in this work are 
part of the spectroscopic Main Galaxy Sample (Strauss 
et al. 2002). 

Rcdshifts (z), de Vaucouleurs model photometry 
{ugriz apparent magnitudes, Fukugita et al. 1996), 
Galactic extinctions, and Pctrosian radii enclosing 50% 
and 90% of the total galaxy light (i?5o and Rqq, respec- 
tively) are taken from the NYU Value-Added Galaxy 
Catalog (Blanton et al. 2005) version of SDSS Data Re- 
lease Four (Adelman-McCarthy et al. 2006). These are 
supplemented with measurements from the SDSS catalog 
archive server* of velocity dispersions {a fib, as measured 
in the SDSS fiber aperture), r-band de Vaucouleurs radii 
(Rdev) and de Vaucouleurs profile axis ratios ((a/b)(iev)i 
as well as the likelihoods of de Vaucouleurs and exponen- 
tial fits to the galaxy light profiles. The catalog archive 
server also provides separate magnitudes of the de Vau- 
couleur and exponential components of each galaxy's 
light (essentially a bulge-disk decomposition) , along with 
the relative contribution of each, which can be used to 
construct model bulge -|- disk magnitudes and colors for 
each galaxy. 

The algorithm for measuring velocity dispersions has 
been updated in Data Release Six (DR6) to correct for 
small systematic biases present in the pre-DR6 mea- 
surements (Bernardi 2007). They are computed using 
the "direct-fitting" method (Burbidge et al. 1961; Rix 
& White 1992), based on fitting stellar templates con- 
volved with gaussian broadening functions. Only a val- 
ues above the ^ 70 km s^^ instnmiental resolution are 
considered reliable. We therefore only use galaxies with 
logo- > 1.86. Velocity dispersion measurement errors are 
typically < 10% (0.04 dcx), based on repeat measure- 
ments of individual SDSS galaxies^. 

^ http;//cas. sdss.org/dr6/en/ 

^ http: / / www.sdss.org/dr6 / algorithms / veldisp.html 
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There is a known problem with SDSS photometry of 
bright galaxies caused by scattered light that biases the 
sky background high. Because of this, the total luminosi- 
ties and radii of bright galaxies are underestimated by the 
SDSS pipeline, as noted by Lauer et al. (2007), Bernardi 
et al. (2007), and Lisker et al. (2007). To correct for 
this effect, we use the results of the simulations reported 
on the SDSS website^ to compute magnitude and radius 
corrections as a function of apparent magnitude, which 
we apply to the galaxies in our sample. We additionally 
limit our sample to galaxies with 15 < r < 18 to avoid 
the worst of the photometry errors caused by this ef- 
fect. The Galactic extinction corrections from the NYU 
Value-Added Catalog are applied to all photometry. 

The galaxy sample presented here is limited to a mod- 
est range in redshift (0.04 < z < 0.08). This has a 
number of advantages. The SDSS Main Galaxy Sam- 
ple is limited to apparent magnitude r < 17.77. For 
0.04 < z < 0.08, our galaxy sample is complete down to 
absolute magnitude Mr = —20.12, while the lower end of 
our redshift range provides data down to Mr > —19. The 
relatively small range in redshifts means that the range 
in lookback times over the entire sample is only At ~ 0.5 
Gyr, which essentially eliminates the effects of evolution 
within our sample. At these redshifts, the SDSS spec- 
tral fibers (1.5" radius) cover a significant portion of the 
galaxy, typically ^ Rdev/2. These spectra are therefore 
not nuclear spectra but sample a substantial fraction of 
the galaxy light. 
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Fig. 1. — Sample galaxies in the color-magnitude diagram. Points 
show the galaxies that meet the sample selection criteria described 
in §2, while contours show all galaxies at 0.04 < z < 0.08. Galaxies 
in the sample are selected to be free of emission lines and to have 
light profiles characteristic of early type galaxies. No explicit color 
selection is applied, yet the resulting galaxy sample consists almost 
exclusively of red sequence galaxies. The remaining color outliers 
are not included in the stacked spectra. See §4.1 for details. 

Throughout the rest of this paper, we will use velocity 
dispersion measurements that have been corrected to 1 /8 
effective radius. Following Bernardi et al. (2003a) and 
J0rgensen et al. (1995), we use a = cr/i6(»'/i(,/|ro)"°'*, 
where To is the circular galaxy radius in arc seconds, 
computed as Tq = RdeV \/ (ci / b) deV , and rfn, is the spec- 
tral fiber radius, 1.5" for SDSS spectra. This allows us to 

® http:/ /www. sdss.org/dr6/start/aboutdr6. html 



provide consistent comparisons with past stellar popula- 
tion work, which typically uses central velocity dispersion 
measurements. In practice, velocity dispersion gradients 
are shallow, and this correction is small (of order ^ 6%). 
The r-band and g-hand magnitudes used throughout this 
paper have been K-corrected to z = 0.0 using the IDL 
code kcorrect vi.lA (Blanton et al. 2003) and converted 
to absolute magnitudes {Mr and Alg) assuming a stan- 
dard ACDM cosmology with A = 0.7, Qm = 0.3 and 
h = 0.7. 

Spectra of the galaxies in our final sample were down- 
loaded from the SDSS DR6 (Adelman-McCarthy 2007) 
version of the SDSS data archive server^. 

2.2. Selection Criteria 

To create a sample of passively evolving galaxies, we 
include only those galaxies with no emission from ionized 
gas. Specifically, we require that the measured equiva- 
lent widths of Ha and [Oii]A3727 emission be below a 2a 
detection threshold, as determined by Yan et al. (2006). 
In this sample, the typical 2a detection threshold corre- 
sponds to Ha and [Oil] equivalent widths of ^ 0.3 A and 
~ 1.7 A, respectively. This criterion excludes galaxies 
with significant ongoing star formation within the cen- 
tral portion of the galaxy covered by the spectral fiber, as 
well as active galactic nuclei (AGN) with optical emission 
hues (Yan et al. 2006). 

Because the SDSS spectra only sample galaxy light 
with the 3" spectral fiber radius, selecting quiescent 
galaxies in this way does not exclude disk galaxies with 
quiescent central bulge populations which harbor star 
formation in their outer regions. These galaxies are not 
truly quiescent galaxies, although they may have quies- 
cent bulges. To limit contamination by these galaxies, 
we follow Bernardi et al. (2003a) in excluding galaxies 
with Rqo/R^o < 2.5 in the i-band (eliminating galax- 
ies without centrally-concentrated light profiles) and for 
which the likelihood of a de Vaucouleurs light profile fit 
is less than 1.03 times the likelihood of an exponential 
fit. These further cuts reduce the sample size by < 10%, 
leaving a total sample of 16,008 galaxies. 

In summary, the galaxies in this study have been se- 
lected from the SDSS Main Galaxy Sample using the 
following criteria: 

1. Redshift range: 0.04 < z < 0.08 

2. Apparent magnitude range: 15 < r < 18 

3. No detected emission in Ha and [Oii]A3727 

4. Concentrated hght profiles: Rqo/R^q > 2.5 in the 
i-band 

5. Likelihood of a de Vaucouleurs light profile > 1.03 
times larger than the likelihoood of an exponential 
light profile 

Throughout the rest of this paper, we refer to the 
galaxies in our sample as quiescent or passive galaxies, 
with the understanding that this categorization refers to 
galaxies which are not actively forming stars. 

In Figure 1, the galaxies used in this analysis are shown 
with respect to the color-magnitude distribution of all 

http://das.sdss.org/DR6-cgi-bin/DAS 
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Fig. 2. — Thumbnail images of a random selection of galaxies from the sample. Circles indicate the size of the 3" SDSS spectral fibers. 
The sample is dominated by early type galaxies, but contains some bulge-dominated spirals as well at a range of orientations. Even where 
a disk is present, the spectral fiber predominantly samples light from the bulge. 



galaxies in this rcdshift range. The gray points show 
the galaxies included in the sample, while the contour 
lines show the underlying color-magnitude distribution 
of all galaxies with 0.04 < z < 0.08 from the SDSS Main 
Galaxy spectroscopic sample. In choosing the sample 
presented here, no explicit color selection was applied, 
yet the resulting sample of galaxies populates a narrow 
red sequence relation between galaxy Mr and color. Fig- 
ure 2 shows thumbnail images for a randomly-selected set 
of galaxies from our final sample. The galaxies all have 
smooth morphologies. About 20% of the galaxies have 



significant disks, although all show evidence for promi- 
nant bulges. In the appendix, we show that the use of 
magnitudes and colors from de Vaucouleur fits to the 
galaxy light profiles (even in galaxies which have a disk 
component) has a negligible effect on our results. 

Of our sample galaxies, 85% have bulge fractions above 
0.8 while only 3% have bulge fractions below 0.5, mak- 
ing the sample galaxies highly bulge-dominated. Figure 
2 shows that the sample contains a modest fraction of 
face-on and edge-on galaxy disks, however these do not 
dominate the light of the galaxies. Furthermore, the disk 
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contribution to the measured spectrum will be substan- 
tially smaller than its contribution to the total galaxy 
light. The 3" SDSS spectral fibers (shown as circles in 
Figure 2) sample roughly 0.3-0.6 of the typical galaxy 
effective radius at these redshifts. In face-on disks, the 
SDSS 3" spectral fiber predominantly samples the galaxy 
bulge, while the spectra of inclined disks will contain 
a slightly increased light contribution from the galaxy 
disks. Edge-on disks in the sample must have negligible 
ongoing star formation in order to satisfy the stringent 
emission line cuts. 

3. THE COLOR-MAGNITUDE RELATION 

The selection criteria defined above result in a sample 
of quiescent galaxies which populate the red sequence 
in the color-magnitude diagram. These galaxies show a 
clear color-magnitude relation (Figure 3, top panel), such 
that more luminous galaxies have redder colors. Data 
points are color-coded by bins in a as labelled in the 
lower figure panels. In the lower panels, the total galaxy 
sample is shown as a cloud of gray data points, with the 
total color-magnitude relation overplottcd as the gray 
dotted line. In each panel, contours indicate the color- 
magnitude distribution of galaxies within a single bin in 
cr, as labeled. These a bins are used throughout this pa- 
per. The vertical dashed line in Figure 3 indicates the 
completeness limit of the SDSS Main Galaxy Sample in 
the redshift range used for this analysis. The solid line at 
Mj. = —19.7 shows where we expect the sample is miss- 
ing more than 50% of the galaxies at these magnitudes. 
Thus the lack of galaxies at fainter magnitudes is likely 
a selection effect rather than a genuine lack of galaxies 
at Mr > -19.7. 

What is notable about the contours in Figure 3 is 
that they show no color-magnitude relation at fixed a, 
in agreement with the work of Bernard! et al. (2005). In 
each panel, which illustrates a narrow range in c, the 
contours in the color-magnitude diagram are horizontal. 
However, galaxies with larger a (red points) are gener- 
ally more luminous and have redder colors than those 
with smaller a (blue and purple points), so that super- 
posing the various a bins results in a net color-magnitude 
relation. 

In contrast, the color-cr relation is readily apparent at 
fixed Mr (Figure 4, top). Here, colored data points in- 
dicate bins in luminosity, with dark green, light green, 
yellow, orange, pink, and purple data points correspond- 
ing to -18.8 > Mr > -20.0, -20.0 > Mr > -20.4, 
-20.4 > Mr > -20.8, -20.8 > Mr > -21.2, -21.2 > 
Mr > -21.6, and -21.6 > Mr > -23, respectively. A 
clear color-cr relation is apparent at fixed Mr (e.g., pur- 
ple data points), with roughly the same slope at the total 
color-cr relation. This color-cr correlation is stronger and 
tighter than the color-magnitude relation shown in Fig- 
ure 3. 

Finally, the bottom panel of Figure 4 shows the a- 

magnitude or Faber- Jackson relation, color-coded by 
galaxy g — r color. Dark blue, medium blue, light blue, 
pale pink, medium pink, and dark pink data points repre- 
sent g-r < 0.74, 0.74 < g-r < 0.76, 0.76 < g-r < 0.78, 
0.78 < g-r < 0.80, 0.80 < g-r < 0.82, and 0.82 < g-r, 
respectively. The cr-magnitude relation is tight for the 
most luminous galaxies, but shows considerable spread 
for faint galaxies. The spread is correlated with galaxy 



color, such that bluer galaxies tend to have lower a at 
fixed luminosity. This may indicate that fainter, bluer 
galaxies have a larger degree of rotational support, re- 
sulting in lower values of a for a given M* . 

It is clear that cr, luminosity, and color are all cor- 
related in early type galaxies. Bernardi et al. (2005) 
argue that the fundamental relations are cr-color and cr- 
magnitude and that the color-magnitude relation results 
from a combination of these two correlations. The color- 
cr relation (Figure 4, top) exists independent of galaxy 
luminosity (i.e., at fixed Mr). Likewise, the cr-magnitude 
relation (Figure 4, bottom) exists independent of galaxy 
color. However, the color-magnitude relation does not 
exist at fixed a (Figure 3), indicating that the cr depen- 
dence of both color and magnitude is what leads to the 
observed color-magnitude relation. 

From Figure 3, it is also clear that galaxies with very 
similar values of cr show a substantial spread in M,.. The 
low cr bins suffer from incompleteness effects, but in the 
higher a bins (e.g., those with log cr > 2.18), it is clear 
that galaxies with a range of only 0.09 dex in log cr span 
roughly 2 mag (0.8 dex) in Mr- Similarly, the substantial 
overlap of the a bins in the color-magnitude diagram 
means that galaxies at fixed Mr have a range of cr. 

These arguments illustrate that optical luminosity L is 
not a good proxy for cr. This is expected from the fact 
that galaxies populate the FP and the Faber- Jackson re- 
lation is not an edge-on projection of the FP (Dressier 
et al. 1987). The two parameters are, however, often used 
interchangeably in the literature as measures of galaxy 
size or mass. We will see in §5.1 that using L versus using 
cr as the controlling variable in stellar population studies 
can lead to substantially different stellar population scal- 
ing relations, because stellar population properties vary 
systematically with L at fixed cr. 

4. DISSECTING THE COLOR-MAGNITUDE RELATION 

Accurate stellar population parameters can only be de- 
rived from high S/N spectra {S/N > 100 A^i; Cardiel 
et al. 1998). As individual SDSS galaxy spectra in our 
sample typically have S/N » 20 A~^, dozens or more 
spectra must be combined to achieve the necessary S/N. 
If galaxies are combined randomly, the average spectra 
would all look very similar. The challenge in this method 
is to identify useful ways of sorting the galaxies in ad- 
vance, in order to do the best possible job of populating 
the full range of galaxy star formation histories, without 
averaging out all the interesting variations. 

In this section, we define a three-dimensional param- 
eter space of cr (which is not sensitive to stellar popula- 
tion effects) along with L and color (which are) and sort 
galaxies into bins in this space. We then stack the spectra 
of galaxies within each bin and measure absorption line 
strengths. In §5, we combine the observed line strengths 
with stellar population models to obtain the fundamen- 
tal stellar population parameters of galaxies along and 
across the color-magnitude relation. 

4.1. Binning by a . L. and color 

In §3, we divided our sample of quiescent galaxies into 
six bins based on cr. To study how stellar populations 
vary at fixed cr, we further divide those same cr bins into 
three bins in color by three bins in L for a total of 54 bins. 
Figure 5 illustrates this process for the logcr = 2.18-2.27 
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Fig. 3. — Top panel: The color-magnitude relation for all galaxies in our sample. Data points are color-coded by a, with purple (red) 
indicating the lowest (highest) a bin. The dashed vertical lines show the completeness limit of the sample, while the solid vertical lines 
show the magnitude at which the sample is missing 50% of the galaxies (i.e., where the sample becomes severely incomplete). Lower panels: 
The color— magnitude relation for galaxies at fixed cr. Gray points show the color-magnitude relation for the total galaxy sample, while 
over-plotted contours show the color-magnitude distribution of galaxies within a limited range in cr, as labeled. Contour lines indicate a 
number density of 5, 10, 20, 40, 60, 80, 100, and 120 galaxies per bin, where bins in Mr and g — r are defined by the tickmarks on the axes. 
The gray dotted line indicates the color-magnitude relation for the total sample. Galaxies with larger values of cr tend to be brighter and 
redder than those at lower tr. However, galaxies at fixed cr span a substantial range in both color and magnitude, so that different bins in 
cr show a large amount of overlap in both color and magnitude. At fixed cr, there is no color-magnitude relation. Combining galaxies from 
all cr bins results in an overall color-magnitude relation. 



bin. Cuts in color and L are defined with fixed width, 
such that each L bin is 0.8 mag wide in Mr and each 
color bin is 0.032 mag wide in g — r. The 3x3 grid of 
sub-bins is centered on the median Mr and g — r for all 
the galaxies in the given a bin. The median values of cr, 
Mr, and g — r for each of the 54 bins are listed in Table 
CI, along with the total number of galaxies in each bin. 
The central sub-bin in color-L space contains more 



galaxies than the outer bins, but all bins are sufficiently 
well populated to produce high S/N stacked spectra. 
Defining the color-L grid in this way excludes color and L 
outliers from each a bin, which ensures that the stacked 
spectra will not be biased by very deviant or misidenti- 
fied objects. 

In the lowest a bins, the magnitude limit of the sample 
presented here means that the fainter galaxies at that a 
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Fig. 4. — Top panel: The color-cr relation, with data points color-coded by galaxy luminosity. Dark green (purple) points show the 
faintest (brightest) luminosity bins. A color-o- relation exists at fixed Mr (e.g., purple data points alone) with roughly the same slope as 
the total color-(T relation. The color-c trend is stronger and narrower than the color-magnitude relation in Figure 3. Bottom panel: The 
cr-magnitude relation, color-coded by galaxy g — r color. Dark blue (dark pink) points show the bluest (reddest) galaxies. The cr-magnitude 
relation is broader for fainter galaxies. This is correlated with galaxy color, such that bluer galaxies tend to have lower a at fixed Mr- This 
may be due to a greater degree of rotational support in fainter and bluer early type galaxies, showing lower <j for the same M*. 



are under-represented. Under the assumption that the 
missing faint galaxies at z ^ 0.08 are not substantiaUy 
different from those at z ~ 0.04 which are included in 
the sample, the fainter bins are not strongly biased but 
are less populated than they ought to be. Comparisons 
between galaxies at different L in the same a bin should 
therefore be robust to completeness effects. However, 
because these a bins are missing the faintest galaxies, 
the median Mr used to locate the color-i grid may be 
biased brighter than the underlying galaxy population at 
that (T. 

4.2. Constructing the stacked spectra 

In each of the 54 bins in <T-L-color space, we combine 
spectra of the galaxies in the bin in order to produce a 
very high S/N average spectrum. Regions around the 
bright skylines at 5577A, 6300A, and 6363A are masked 
in the individual galaxy spectra, which are then shifted to 



the restframe. The individual spectra are normalized so 
that all galaxies in the bin contribute equally to the aver- 
aged spectrum. The normalization uses the median flux 
in the 4100-5000A range without modifying the contin- 
uum shape, as the flux-calibrated spectral shape affects 
the Lick index measurements. In practice, because the 
galaxies in a given bin have very similar colors, the shape 
of the continuum is nearly identical for all galaxies in the 
bin. 

The stellar absorption features in higher a galaxies 
are effectively observed at lower resolution than those in 
lower (T galaxies, due to the increased intrinsic Doppler 
smoothing of the lines from stellar motions within the 
galaxy. This intrinsic smoothing has non-negligible ef- 
fects on the measured absorption line strengths. In order 
to compare all galaxies at the same effective resolution, 
we smooth lower-u galaxies up to match the highest-a 
galaxies in our sample, a « 300 km s~^. The algorithm 
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Fig. 5. — Method used to bin galaxies by L and color at fixed a. 
The color-magnitude relation for the total quiescent galaxy sample 
in this typical a bin is indicated by the gray contours. Overplottcd 
as gray points are all galaxies with 2.18 < log cr < 2.27. Within 
this (T bin, there is no color-magnitude relation. The median point 
in Mr and g — r is shown as the black cross in the center. Galaxies 
in this (7 bin are divided into a 3x3 grid in color-magnitude space, 
with the grid centered on the median point and with fixed grid 
spacing of 0.8 mag in Mr and 0.032 mag in g — r. Defining the 
grid in this way eliminates color and L outliers which may bias the 
average properties of the bins. Galaxies in each of the six cr bins 
are sorted in this way into nine sub-bins in color-L space, and the 
galaxies in each sub-bin are stacked as described in §4.2 to produce 
high S/N mean spectra for all 54 bins in cr-L-color space. 



we use to stack the galaxy spectra first averages the un- 
smoothed spectra of galaxies within narrow a sub-bins, 
with a routine that rejects very deviant pixels. The av- 
erage clean sub-bin spectra are then smoothed to the 
maximum a of the sample 300 km s^^) and coadded. 
The result is a stacked spectrum for each cr-L-color bin 
from which outlier pixels have been rejected and which 
has been smoothed to an effective resolution identical to 
all the other bins. The error spectrum is also computed 
from the individual SDSS galaxy error spectra by adding 
the individual pixel errors in quadrature to produce an 
error spectrum for each stacked spectrum. 

4.3. Lick Index absorption strengths 

We measure the full set of Lick indices in each of the 
stacked spectra. These include the Balmer lines H/3, 
a^F, and HSp (also broad definitions of Hja and USa), 
and a set of Fe-dominated hues (Fe4383, Fe4531, Fe5015, 
Fe5270, Fe5335, Fe5406, Fe5709, and Fe5782), as weU 
as numerous lines that are sensitive to abundances of 
elements other than Fe (Mgi, Mg2, Mg b, CNi, CN2, 
C24668, Ca4227, Ca4455, NaD, TiOi, and TiOa). The 
index definitions are taken from Worthey et al. (1994) 
and Worthey & Ottaviani (1997), and line strengths are 
measured using the Lick_EW code that is available on- 
line^ as part of the EZ_Ages code package (Graves & Schi- 
avon 2008). The Lick_EW code also computes statistical 
errors for each Lick index based on the error spectra, us- 
ing equations 33 and 37 of Cardiel et al. (1998). Because 
the Lick indices are defined as equivalent widths, they 
are insensitive to dust extinction effects. 

® http:/ /www. ucolick.org/~gravcs/EZ_Ages. html 



All the stacked spectra have been smoothed to a = 300 
km s"^. Combined with the intrinsic resolution of the 
SDSS spectrograph, the spectra are at lower resolution 
than the Lick/IDS resolution at which the Lick indices 
are defined. A correction must be applied to the mea- 
sured line strengths to bring them back onto the Lick 
index system. These corrections are included in the 
Lick_EW computation, based upon smooth single stel- 
lar population spectra, as given in Schiavon (2007) Table 
A2a. All of the spectra have been smoothed to the same 
a and corrected in the same way so that uncertainties in 
the velocity dispersion correction should not affect rela- 
tive lines strength measurements. The corrected Lick in- 
dex measurements and statistical errors for each of the 54 
stacked spectra are given in Table C2. We do not attempt 
to match the zeropoints of the SDSS spectra to those of 
the Lick system as defined by Schiavon (2007). This 
should only introduce minor zero-point uncertainties be- 
cause the Schiavon models are based on flux-calibrated 
spectra and the SDSS spectra are also fiux calibrated (see 
Schiavon 2007, section 2.2.2). 

4.4. Stellar population modelling 

We compare the Lick index measurements for each of 
the stacked spectra to the stellar population models of 
Schiavon (2007) to derive fundamental stellar population 
parameters from the line strengths. The Schiavon (2007) 
models are single burst models which include the effect of 
variable abudance ratios by combining theoretical stellar 
isochrones (Girardi et al. 2000) with a library of empirical 
stellar spectra (Jones 1999) and individual line strength 
sensitivities to elemental abundance variations computed 
from theoretical stellar atmospheres (Korn et al. 2005). 
With these models, a set of Lick index measurements can 
be used to determine the mean luminosity-weighted stel- 
lar population age, iron abundance ([Fe/H]), and abun- 
dance ratios for the elements Mg, G, N, and Ca ([Mg/Fe], 
[C/Fe], [N/Fe], and [Ga/Fe]) using the algorithm de- 
scribed in Graves & Schiavon (2008) and implemented 
in the publicly available^ code EZ_Ages. Briefly, the 
Graves & Schiavon (2008) method determines a fiducial 
mean age and [Fe/H] from a combination of Balmer and 
Fe-dominated indices, then uses other indices which are 
sensitive to Mg, C, N, and Ga to adjust the element abun- 
dance ratios until the model index predictions match the 
ensemble of Lick index data. The basic modelling pro- 
cess has not changed since the original analysis of red se- 
quence galaxies with LINER-like emission in Graves et al. 
(2007); subsequent minor updates to the EZ_Ages code 
have improved the robustness of the code at runtime and 
produce results that match those of Graves et al. (2007) 
within the quoted errors. EZ_Ages estimates stastical er- 
rors for each of the stellar population parameters, based 
on the measurements errors in the Lick indices. The er- 
rors in these various parameters are not independent; an 
analysis of the correlated errors is presented in Figure 3 
of Graves & Schiavon (2008). The reader is referred to 
that article for additional details on the age and abun- 
dance fitting process. 

Graves & Schiavon (2008) have shown that this fitting 
method is robust to the choice of Lick indices used in 
the analysis. Here we use the "standard set" of Lick in- 

^ http:/ /www. ucolick.org/~graves/EZ_Agcs. html 
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Fig. 6. — Mean stellar population parameters for each of the 
54 stacked galaxy spectra, shown as a function of median a for 
the galaxies in each bin. The error bars in the lower right cor- 
ner of each panel show the median statistical error in the derived 
population parameters, based on the uncertainties in the Lick in- 
dex measurements. The data have been separated into those with 
small statistical errors (black) and those with larger statistical er- 
rors (gray), as indicated by the corresponding error bars. "Large" 
statistical errors are defined as those more than two standard de- 
viations larger than the median error. Mean luminosity-weighted 
stellar age, [Fe/H], [Mg/H], and [Mg/Fe] all increase with increas- 
ing (T. Dashed lines show linear least squares fits of the stellar 
population properties onto log ct, computed using the black data 
points only. The Spearman rank correlation coeffecients (r^) of 
the correlations are specified in each panel. The scatter in mean 
age and [Fe/H] is significantly larger than that expected due to 
measurement errors, implying genuine variation in these popula- 
tion parameters at fixed u. The spread in age is small for the 
highest-(T galaxies, all of which show old ages, while the lowest-cr 
galaxies cover a substantial factor in mean age. There is only mod- 
est scatter in either [Mg/Fe] or [Mg/H] at fixed a. The correlation 
bewteen cr and [Mg/H] is particulary strong, as indicated by the 
very high value of r^. 



dices from Graves & Schiavon (2008) — H/3, (Fe) (which 
is an average of the Fe5270 and Fe5335 indices), Mg &, 
C24668, CNi, and Ca4227~with the EZ_Ages code to 
fit a mean luminosity- weighted age, [Fe/H], [Mg/Fe], 
[C/Fe], [N/Fe], and [Ca/Fe] for each of the 54 average 
galaxy spectra. We can thus track each of these stellar 
population parameters as a function of cr, L, and color 
for the galaxies in our sample. 

In this analysis, we focus on the stellar population 
age, [Fe/H], [JVIg/H], and [Mg/Fe], leaving the analysis 
of other abundance ratios for future work. The results 
of the stellar population analysis for the stacked spec- 
tra are presented in Table C2. Of the a-elements, Mg 
is the only one that is relatively straightfoward to mea- 
sure in optical galaxy spectra. For the purposes of this 
type of stellar population work, it is typically assumed 
that all a-elements scale together and thus that [Mg/Fe] 
is equivalent to [a/Fe]. There is some observational evi- 
dence that a elements may not always vary in lock-step 
(e.g., Fulbright et al. 2007; Humphrey & Buote 2006) 
but, lacking an alternative, we will use [Mg/Fe] as an es- 
timate of the total [a/Fej. The Schiavon (2007) models 
are computed at fixed [Fe/H] rather than at fixed total 
metallicity ([Z/H]), so that enhancements in individual 
elements also change the total metallicity of the galaxy. 

It is important to keep in mind that the stellar pop- 



FlG. 7. — Mean stellar population parameters as a function of 
galaxy L. The results of the stellar population analysis are plotted 
against the median Mr for the galaxies in each bin. Statistical 
errors in the derived population parameters are shown in the lower 
right corner of each panel. The data have been separated into 
those with small statistical errors (black) and those with larger 
statistical errors (gray), as indicated by the corresponding error 
bars. The Spearman rank correlation coefficient [ts) is given for 
each relation. The strong increases in mean age and [Mg/Fe] visible 
in Figure 6 when plotted against a are considerably washed out 
when plotted against Afr, while the trend between [Fe/H] and Mr 
is much tighter than the (T-[Fe/H] relation. The [Mg/H]-M,- is 
still strong, although Ts indicates it is slightly weaker than the 
(T-[Mg/H] correlation. 

ulation parameters based on modelling absorption line 
indices give ages and abundances averaged over the en- 
semble of stars in the galaxy. The contribution of each 
star to the average age or abundance measurement is 
proportional to the light contributed at the wavelength 
of the absorption feature in question, and also to the 
equivalent width of the particular absorption feature of 
that star. When comparing ages of galaxies, a younger 
mean luminosity-weighted age in one galaxy compared 
to another does not necessarily mean that all of the stars 
in the galaxy are younger: a galaxy may contain a frac- 
tion of young stars that skew the average age to lower 
values. The luminosity-weighted ages presented here are 
also averaged over the ensemble of galaxies in the bin 
and should thus be treated as a statistical description of 
galaxies in that bin rather than as "true" ages for all of 
the stars in an individual galaxy. 

5. STELLAR POPULATIONS ON THE RED SEQUENCE 

5.1. Age, [Fe/H], [Mg/H], and [Mg/Fe] as Functions of 
a, L, and Color 

Figure 6 shows the age, [Fe/H], [Mg/H], and [Mg/Fe] 
results of the stellar population analysis as a function 
of galaxy cr. The stellar population parameters derived 
from the stacked spectra are plotted against the median 
value of cr for the ensemble of galaxies in each of the 54 
bins. The median statistical errors in the derived pop- 
ulation parameters are shown in the lower right corner 
of each panel. The Spearman rank correlation coeffi- 
cient (rg)^^ is shown for each pair of parameters, in- 
dicating the strength of the correlation. Age, [Fe/H], 

For a given pair of parameters, rs indicates the extent to 
which they produce the same rank ordering of objects, with in- 
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Fig. 8. — Mean stellar population parameters as a function of 
galaxy color. The results of the stellar population analysis are 
plotted against the median g — r color for the galaxies in each bin. 
Statistical errors in the derived population parameters are shown in 
the lower right corner of each panel. The data have been separated 
into those with small statistical errors (black) and those with larger 
statistical errors (gray), as indicated by the corresponding error 
bars. The Spearman rank correlation coefhcienct (rs) is given for 
each relation. Here, mean age increases with color, similar to the 
(T— age trend observed in Figure 6 although with somewhat smaller 
scatter in log age at fixed g — r than was seen at fixed a. However, 
[Fe/H] shows only a very weak trend with color, while [Mg/Fe] 
and [Mg/H] increase significantly in redder galaxies. Neither the 
[Mg/Fe]— color nor [Mg/H]-color trends are as strong as the those 
with (T, as indicated by r^. 



[Mg/H], and [Mg/Fe] all show correlations with cr, such 
that higher-cr galaxies have older mean ages, higher iron 
abundances, higher magnesium abundances, and higher 
levels of Mg-enhancement, which strongly suggests that 
higher-(T galaxies also have higher total metallicities. 
These results are in agreement with numerous previous 
authors (e.g., Bernardi et al. 2003b; Thomas et al. 2005; 
Nelan et al. 2005; Smith et al. 2007a; see summary in 
Table 8 of Graves et al. 2007), although these trends are 
not universally observed in all early type galaxy samples 
(e.g.. Kelson et al. 2006; Trager et al. 2008). 

The [Mg/H]-cr and [Mg/Fe]-cr relations are particularly 
tight compared to the typical measurement errors. This 
implies that at fixed a [Mg/H] and [Mg/Fe] vary only 
weakly with the other two global parameters studied here 
(L and color). At high cr, [Mg/H] shows almost no varia- 
tion at fixed a. A linear least-squares fit of [Mg/Fe] onto 
a gives a slope of 0.34, which is consistent with previous 
results. 

Unlike [Mg/H] and [Mg/Fe], both age and [Fe/H] show 
scatter at fixed a that is large compared to the statistical 
error bars. For [Fe/H], the scatter is roughly the same 
at all (T, whereas the spread in mean ages at fixed a ap- 
pears to be larger for lower-cr galaxies. It is noteworthy 
that some 1ow-ct galaxies have mean ages nearly as old 
as some of the highest-cr galaxy bins. Thus any model 
for galaxy formation that predicts a correlation between 

dicating no correlation and 1 indicating perfect correlation. The 
sign of indicates whether the correlation is positive or negative. 
The definition of Ts makes no assumptions about the functional 
form of the correlation and allows robust comparison of correlation 
strengths between different correlations. 



galaxy a and galaxy age must also be able to accommo- 
date low-CT galaxies with old ages. 

Because a correlates with both L and color (see §3), 
one might expect to see similar correlations between stel- 
lar population properties and L or color as were seen for 
cr. Figure 7 shows the derived stellar population param- 
eters as functions of L, while Figure 8 shows the same as 
functions of galaxy color. Interestingly, when the stellar 
population properties are plotted against L instead of cr, 
the correlations with age and [Mg/Fe] are nearly erased, 
while the correlation with [Fe/H] is much tighter. Con- 
versely, when they are plotted against color, the corre- 
lation with mean stellar age is slightly tighter than with 
cr, while the [Fe/H] and [Mg/Fc] correlations are weaker 
than with a. Of the four stellar population parameters, 
[Mg/H] alone shows relatively strong correlations with 
all three global parameters, although it is most strongly 
correlated with a. 

Figures 6-8 show that, although cr, L, and color are 
all correlated with one another, the various stellar pop- 
ulation parameters show significantly different behav- 
ior with each of the global parameters. Residuals from 
the cr-stellar population relations must correlate with L 
and color in such a way as to erase some trends and 
strengthen others when L or color are used instead of a. 
This correlation of residuals immediately suggests that 
quiescent galaxies populate a multi-dimensional stellar 
population parameter space. To explore these effects in 
detail, we must investigate the dependence of stellar age, 
[Fe/H], [Mg/H], and [Mg/Fe] on aU three global param- 
eters simultaneously. 

5.2. Mapping Stellar Populations Across the Red 
Sequence 

Figure 9 shows contours of stellar population age as 
a function of ct, L, and color. The six panels contain 
data for our six standard a slices corresponding to the 
panels in Figure 3. Black and gray data points in each 
panel show the color- magnitude relation at that cr, with 
gray points indicating galaxies in the cr range which were 
excluded from the stacked spectra. 

Plotted over the color-magnitude relations in each bin 
are color contours representing the mean luminosity- 
weighted ages derived from the stacked spectra. The age 
contours were constructed as follows: age values from 
the stacked spectra in each a range were plotted at the 
median values of and g — r for each bin, forming a 
3x3 grid of age values in color-magnitude space. Val- 
ues of age were linearly interpolated (in log age) between 
the nine grid points to produce age contours across the 
color-magnitude relation. The color bar in the lower right 
panel indicates the scaling of the age contours, with the 
lowest and highest values indicated. The dashed lines at 
Mr = —21 and g — r — 0.75 are the same in all panels 
and exist merely to guide the eye. As cr increases, galax- 
ies become systematically both brighter and redder with 
respect to the dashed lines. 

The age contours in Figure 9 illustrate the same trends 
that were visible in Figure 6. As cr increases from low 
(upper left panel) to high values (lower right panel), the 
average age indicated by the contours increases. How- 
ever, in each slice in a, a range of ages exists. The spread 
in ages at fixed cr is large in the 1ow-ct bins and becomes 
increasingly smaller toward the highest-cr bins. 
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Fig. 9. — Mean luminosity-weighted age as a function of a, L, and color. Each panel shows the color-magnitude diagram for a slice in 
cr, corresponding to the six a bins in Figure 3. Black points represent galaxies that fall within 1.2 mag of the median and within 0.048 
mag of the median g — r for the a bin. Gray points show galaxies that fall outside these boundaries and which are excluded from the 
stacked spectra, as described in §4.1. Overplotted are contours of mean luminosity-weighted galaxy age, derived from the stellar population 
modelling in §4.4. The dashed lines are reference lines to guide the eye when comparing between a bins. Stellar population age increases 
with increasing cr, but shows substantial scatter within a fixed cr slice. There is more age variation amongst low-cr galaxies than in high-cr 
galaxies. At fixed cr, lines of constant age run roughly diagonally across the color-magnitude diagram, with younger galaxies having bluer 
colors and brighter luminosities, while older galaxies have redder colors and fainter luminosities. When the various a slices are stacked 
together into the observed red sequence, the decrease in age with L at fixed cr counter-acts the overall increase of age with cr, erasing the 
age trend when L is used to parameterize galaxies instead of cr (as in Figure 7). See §5.2 for details. 



The contours also illustrate the behavior of stellar pop- 
ulation age at fixed a across the color-magnitude rela- 
tion. Lines of constant age thus run roughly diagonally 
in color-magnitude space at fixed tr, with young ages in 
the lower right corner (bright and blue) and old ages in 
the upper left corner (faint and red). This trend is ap- 
proximately repeated in each a slice, although at high 
a the dynamic range is less. The variation at fixed a 
cannot be explained by errors in a that result in galaxies 
being assigned to the wrong stacked spectrum because 
this would result in the opposite effect: galaxy contami- 
nation from higher-(T bins would tend to make the more 
luminous galaxies appear older, not younger. 

Figure 10 shows the same slices through the color- 
magnitude relation, this time with color contours indi- 
cating [Fe/H]. Again, the trends of Figure 6 are visible: 
[Fe/H] increases with a from the lowest-u bins to the 
highest, and a range in [Fe/H] exists within each slice 



in a. However, the variation in [Fe/H] at fixed a is 
markedly different from the age variation seen in Fig- 
ure 9. Lines of constant [Fe/H] run nearly vertically in 
the color-magnitude diagram for fixed cr, such that more 
luminous galaxies have higher [Fe/H]. 

The contours of [Mg/H] shown in Figure 11 are similar 
to those in [Fe/H] for the same cr slices through the color- 
magnitude relation. In addition to the total increase in 
[Mg/H] with cr, the brightest galaxies at fixed cr have the 
highest values of [Mg/H]. However, [Mg/H] shows less 
variation at fixed a than does [Fe/H], particularly in the 
high-cr bins. This is as expected from Figure 6. Although 
the spread in [Mg/H] at fixed cr seen here and in Figure 6 
is small. Figure 1 1 strongly suggests that these variations 
are in fact real because they are consistent across the cr 
bins. 

Finally, Figure 12 shows cr slices overlaid with contours 
of constant [Mg/Fe], directly comparing the differences 
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Fig. 10. — [Fe/H] as a function of a, L, and color. As in Figure 9, panels show the color-magnitude diagram for slices of fixed a. 
Black points show galaxies included in the stacked spectra, while gray points show color and M,- outliers that have been excluded from 
this analysis. Overplotted are contours of [Fe/H] derived from the stellar population modelling of §4.4. [Fe/H] increases with increasing 
a. Lines of constant [Fe/H] run approximately vertically in color-magnitude space at fixed tr, with the most luminous galaxies being the 
most Fe-rich. When the various a slices are stacked together into the observed color-magnitude relation, the increase in [Fe/H] with L 
at fixed cr reinforces the increase in [Fe/H] with a, causing the L-[Fe/H] relation to be tighter and steeper than the (T-[Fe/H] relation. 
Comparisons with Figure 9 show that, at fixed cr, age and [Fe/H] are roughly anti-correlated, helping to preserve the tightness of the total 
color-CT relation. 



in Fe and Mg enrichment. Again, [JVIg/Fe] increases with 
increasing a, as expected from Figure 6. Like [IVIg/H], 
the variation at fixed a is mild and therefore somewhat 
more difficuh to categorize than the variation in age and 
[Fe/H], but broadly [Mg/Fe] behaves similarly to age, 
showing larger enhancements in fainter galaxies at fixed 
a. Because JVIg is produced by supernovae Type II (SNe 
II) on short timescales while Fe is produced by super- 
novae Type la (SNe la) on longer timescales, the ob- 
served overabundance of Mg relative to the solar abun- 
dance pattern is typically interpreted to indicate a short 
timescale for star formation, leaving the stellar popula- 
tion of the galaxy enhanced in SN II products but not in 
SN la products. 

These maps of age, [Fe/H], [IVIg/H], and [Mg/Fe] 
across the color-magnitude relation confirm that qui- 
escent galaxies form a multi-parameter family in stel- 
lar population properties. Age, [Fe/H], [Mg/H], and 
[IMg/Fe] all increase with a, but show variations at fixed 
a depending on their L. Of the four stellar population 



parameters presented in Figures 9-12, only age shows a 
strong dependence on color as well as L. Variations in 
age and [IVIg/Fe] with L at fixed a are similar (both in- 
crease toward fainter galaxies) and are opposite to the 
variations in [Fe/H] and [Mg/H] (both of which increase 
toward brighter galaxies). This suggests that the old- 
est galaxies in the universe formed their stars over short 
timescales, while younger quiescent galaxies experienced 
more extended star formation. 

Using these stellar population maps, the differing 
trends with a and with L illustrated in Figures 6 and 
7 can be understood by superposing the various a slices 
of the color-magnitude relation. Although we have seen 
that age and [Fe/H] both increase with cr. Figures 9 and 
10 show that age and [Fe/H] exhibit opposite behavior 
from one another at fixed cr, namely that age decreases 
and [Fc/H] increases for brighter galaxies at fixed cr. Be- 
cause of this difference, superposing the various a slices 
to form the total color-magnitude relation acts to rein- 
force the cr- [Fe/H] relation into an even tighter L-[Fe/H] 
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Fig. 11. — [Mg/H] as a function of a, L, and color. Panels show the color-magnitude diagram for slices of fixed a. Black points 
shows galaxies included in the stacked spectra, while gray points show color and Mr outliers that have been excluded from this analysis. 
Overplotted are contours of [Mg/H] derived from the stellar population modelling of §4.4. Mg is used to trace a-elements, thus contours 
of [Mg/H] here are a proxy for contours of [«/H]. [Mg/H] increases with increasing cr, showing the least variation at fixed a of any of 
the stellar population parameters. In the lower-cr bins, [Mg/H] behaves similarly to [Fe/H], with lines of constant [Mg/H] running nearly 
vertically in color-magnitude space. 



relation. At the same time, the age trends at fixed a 
counteract the cr-age relation, resulting in almost no L- 
age relation. The explains why stellar population studies 
as a function of L or M^,,^^ find weak or non-existant L- 
age relations (e.g., Kuntschner & Davies 1998; Terlevich 
& Forbes 2002; Gallazzi et al. 2005), while stellar popu- 
lation studies as a function of a turn up significant cr-age 
correlations (e.g., Bernardi et al. 2003b; Thomas et al. 
2005; Nelan et al. 2005; Smith et al. 2007a; Graves et al. 
2007). 

In this scenario, although the L-[Fe/H] correlation is 
stronger than the (T-[Fe/H] correlation, we interpret the 
(T-[Fe/H] trend as the primary relation, with the in- 
creased tightness of the L-[Fe/H] relation being caused 
by correlated residuals from the a-L and (T-[Fe/H] rela- 
tions. In the next section, we use principal components 
analysis to show that the family of stellar populations 
in quiescent galaxies is nearly two-dimensional, with the 
first dimension parameterized by a and the second di- 

M» is closely related to L because red sequence galaxies span 
only a very limited range in Mf,/L 



mension parameterized by correlated residuals from the 
various cr-L, tr-color, and tr-stellar population trends. 

5.3. Principal Components Analysis and Stellar 
Population Residuals 

To quantitatively explore the multi-dimensional space 
of global and stellar population parameters illustrated 
in the previous section, we have performed a principal 
components analysis (PGA; see, e.g., Faber 1973; Trager 
et al. 2000) in the seven dimensional space parameterized 
by cr, L, color, age, [Fe/H], [Mg/H], and [Mg/Fe]. All pa- 
rameters are "standardized" to have a mean of zero and 
a variance of one before performing the PGA. The results 
are shown in Table 1. Although the parameter space is 
nominally seven dimensional, there are only six principal 
components (PGs) because [Mg/H] = [Mg/Fe] + [Fe/H]. 
The first two PGs account for 91% of the variance in the 
population, indicating that the quiescent galaxy popula- 
tion can be well-described by a two-dimensional hyper- 
plane. 

The first PG (PGl) accounts for a fuU 70% of the ob- 
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Fig. 12. — [Mg/Fe] as a function of a, L, and color. Panels show the color-magnitude diagram for slices of fixed a. Black points 
shows galaxies included in the stacked spectra, while gray points show color and Air outliers that have been excluded from this analysis. 
Overplotted are contours of [Mg/Fe] derived from the stellar population modelling of §4.4. Mg is used to trace a-elements, thus contours 
of [Mg/Fe] here are a proxy for contours of [a/Fe]. [Mg/Fe] increases with increasing cr. There is limited variation in [Mg/Fe] at fixed cr. 
A comparisons with Figure 9 shows that [Mg/Fe] varies similarly to age, such that fainter galaxies at fixed a are older and have higher 
[Mg/Fe]. This variation runs in the opposite sense of variation in [Fe/H] at fixed cr, which increases with increasing luminosity. When 
the various cr slices are stacked together in the observed red sequence, the decrease of [Mg/Fe] with increasing L at fixed cr counter-acts 
the trend of increasing [Mg/Fe] with increasing cr, thus the cr-[Mg/Fe] relation is substantially weakened when L is used to parameterize 
galaxies instead of cr 



served variation and is comprised of positive contribu- 
tions from all of the input parameters. This iUustrates 
quantitatively that aU the parameters presented here are 
correlated with one another. We saw in §3 that the color- 
magnitude relation is the result of fundamental cr-color 
and (T-magnitude relations. It therefore seems reasonable 
to treat a as the primary global parameter behind PCI. 

The second PC (PC2) encompasses another 21% of the 
variance and includes significant contributions from all 
parameters except for tr. PC2 therefore represents the 
correlated residuals from the mean relations between a 
and the other parameters. The contributions of g—r, age, 
and [Mg/Fe] to PC2 are positive while those of L, [Fe/H], 
and [Mg/H] are negative. This indicates that at fixed a, 
g — r, age, and []VIg/Fe] are correlated with one another 
and are anti-correlated with L, [Fe/H], and [IVIg/H]. Of 
the global parameters, L contributes the most to PC2. 

5.4. Stellar Population Residuals from a Relations 



Residuals in the four stellar population parameters 
(A log age, A [Fe/H], A [Mg/H], and A [Mg/Fe]) are de- 
fined as the difference between the observed values and 
the mean relations with cr, as indicated by the dashed 
lines in Figure 6. The correlated stellar population resid- 
uals associated with PC2 are shown in Figure 13. In each 
panel, high and low S/N data are indicated by filled and 
open circles, respectively. The error ellipses in the lower 
left corner of each panel indicate the stastical errors from 
the stellar population modelling process, with solid black 
(dashed gray) lines showing the median statistical errors 
for high (low) S/N data. Measurement errors in the in- 
dex absorption line strengths produce correlated errors 
in the derived stellar population parameters. These are 
determined from Monte Carlo simulations in Graves & 
Schiavon (2008) and are indicated by the error ellipse 
orientation. In all panels, residuals are color-coded by 
luminosity residuals; the green circles represent the cen- 
tral bin in L at fixed a (centered on the median L for 
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TABLE 1 
Principal Component Analysis 





Parameter 


PCI 


PC2 


PC3 


PC4 


PC5 


PC6 


log c 


0.43 


0.03 


0.31 


0.37 


-0.53 


0.55 


logL^ 


0.37 


-0.40 


0.04 


0.59 


0.58 


-0.14 


g-r 


0.40 


0.23 


-0.62 


-0.25 


0.32 


0.50 


log Age 


0.32 


0.53 


-0.33 


0.32 


-0.27 


-0.57 




Fe/H] 


0.32 


-0.55 


-0.22 


-0.29 


-0.33 


-0.21 




Mg/H] 


0.43 


-0.21 


0.14 


-0.40 


-0.08 


-0.21 




Mg/Fe] 


0.36 


0.39 


0.59 


-0.34 


0.31 


-0.11 


Eigenvalue 


4.88 


1.49 


0.35 


0.18 


0.07 


0.03 


Percentage of variance 


69.7 


21.3 


5.0 


2.6 


1.0 


0.4 


Cumulative percentage 


69.7 


91.0 


96.0 


98.6 


99.6 


100.0 



Note. — Parameters are "standardized'' to have a mean of zero and a 
variance of one before performing PCA. Although seven parameters are re- 
ported here, there are only six principal components because [Mg/H] = [Fe/H] 
+ [Mg/Fe]. 

^PCA is performed using log L rather than to eliminate confusion due to 
the sign convention of the magnitude scale. 



that (T, as in Figure 5) while blue and red circles indicate 
brighter (Ai « —0.8 mag in r) and fainter (Ai » +0.8 
mag in r) luminosity bins. 

Panel (a) shows an anti-correlation between A log age 
and A [Fe/H], such that older galaxies are typically Fc- 
poor compared to younger galaxies at the same a . This 
anti-correlation at fixed a has been reported by numer- 
ous previous authors (e.g., Worthey et al. 1995; Colless 
et al. 1999; J0rgensen 1999; Trager et al. 2000) and was 
quantified by Trager et al. (2000) as the "metallicity hy- 
perplane" . The anti-correlation conspires to keep the 
color-cr relation tight (c.f.. Figure 4). In contrast to the 
observed age-[Fe/H] anti-correlation at fixed a, we saw 
in Figure 7 that at fixed L, there is significant scatter in 
age but almost no scatter in [Fe/H] (although see Pog- 
gianti et al. 2001^^). This implies that age differences 
will produce color variations at fixed L which cannot 
be offset by correlated variations in [Fe/H]. Thus, the 
color-magnitude relation should show more color varia- 
tion at fixed L than is observed in the color-cr relation at 
fixed (7, consistent with Figures 3 and 4. Furthermore, 
the observed color spread in the color-magnitude rela- 
tion should correlate with age, such that bluer galaxies 
are younger than redder galaxies at fixed L. This is con- 
sistent with Figure 9 and also with the results of Cool 
et al. (2006). 

In addition to the age-[Fe/H] anti-correlation evident 
in Figure 13a, a similar anti-correlation is present be- 
tween age and [Mg/H] (Figure 13c), although it is sig- 
nificantly weaker. A further anti-correlation is observed 
between [Fe/H] and [Mg/Fe] (Figure 13b), such that Fe- 
poor galaxies are more Mg-enhanced than their compari- 
tively Fe-rich counterparts at fixed a. Finally, Figure 13d 
shows a positive correlation between age and [Mg/Fe], 

Unlike the results shown in Figure 7, Poggianti et al. (2001) 
find a range of metallicity at fixed L, which they observe to be anti- 
correlated with age. Their results may differ from ours because 
we have averaged over a large number of galaxies in each stacked 
spectrum, washing out genuine [Fe/H] variation at fixed L and 
(T, or the observed age-metallicity correlation in their data may 
be due to an underestimate of the correlated errors in the stellar 
population modelling process. 



such that older galaxies are more Mg-enhanced with re- 
spect to the solar abundance pattern than their younger 
counterparts at fixed a. 

In panels a-c, the slopes of the observed anti- 
correlations are similar to those expected from correlated 
errors in stellar population modelling, as indicated by the 
error ellipses. However, the observed anti-correlations 
cannot be due to correlated statistical errors for sev- 
eral reasons. Firstly, the typical statistical errors are too 
small to produce the observed spread in stellar popula- 
tion properties. More conclusively, if the observed anti- 
correlations were due to measurement errors, they would 
not be systematically correlated with AL. The observed 
L-dependence of stellar population properties at fixed u 
argues strongly that the anti-correlation is real. In ap- 
pendix B, we verify that the observed anti-correlations 
are not due to systematic effects in the stellar population 
modelling process caused by using single burst models to 
fit integrated galaxy spectra, which are almost certainly 
not single burst populations. 

It is evident from Figure 13 that, not only are the 
various stellar population residuals correlated with one 
another, they are also strongly correlated with AL, as 
quantified by the PCA analysis in the previous section. 
Collecting all parameters together, the residual trends 
can be summarized as follows: at fixed a, galaxies fainter 
than the median are older, Fe-poor, somewhat Mg-poor, 
and more Mg-enhanced than typical galaxies at that a, 
while galaxies brighter than the median are younger, Fe- 
rich, somewhat Mg-rich, and less Mg-enhanced. 

6. DISCUSSION 

We have shown that (1) stellar population age, [Fe/H], 
[Mg/H], and [Mg/Fe] all increase with increasing galaxy 

(T (Figure 6), and (2) that the star formation histories of 
galaxies vary systematically at fixed a such that fainter 
galaxies are older, Fe-poor, and Mg-enhanced c;ompared 
to their bright counterparts at the same a (Figure 13). 
The first of these results is in agreement with a substan- 
tial number of earlier works (e.g, Bernardi et al. 2003b; 
Thomas et al. 2005; Nelan et al. 2005; Smith et al. 2007a; 
Graves et al. 2007). In particular, the increase in age with 
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Fig. 13. — Correlated residuals in stellar population parameters at fixed a. Residuals Alog Age, A[Fe/H], A[Mg/H], and A[Mg/Fe] were 
computed by subtracting off the trends with a (dashed lines in Figure 6) from the observed values. Data are color-coded by luminosity bin 
(AL) as indicated. Filled (open) circles show high (low) S/N data, while solid (dashed) error ellipses show the corresponding correlated 
errors from stellar population modelling. The Spearman rank correlation coefficient (ts) is indicated for each pair of residuals. Residuals 
from the mean relations with a are related such that residuals in age and [Mg/Fe] are correlated with one another and are anti-correlated 
with residuals in [Fe/H] and [Mg/H]. These correlated residuals are also correlated with AL, such that the brighter galaxies (blue points) 
at a given a have lower ages, lower [Mg/Fe], higher [Fe/H], and higher [Mg/H] than their fainter counterparts (red points). Although the 
residuals are in many cases correlated in the same direction as the correlated errors from stellar population modelling, the relation between 
the stellar population residuals and AL could not be produced by stellar population modelling errors and indicates that the residual 
correlations are real. 



(7 is consistent with "archeological downsizing" (Thomas 
et al. 2005), in which massive galaxies form their stars 
predominantly at early times, while lower mass galaxies 
show evidence of younger stellar populations. 

The anti-correlation of age and [Fe/H] at fixed a shown 
in this analysis is not a new result and is in agreement 
with the metallicity hyperplane of Trager et al. (2000) 
and with the recent work of Smith et al. (2007b) for in- 
dividual galaxies. However, an important new result of 
this work is that this anti- correlation is correlated with 
galaxy luminosity, such that the more luminous galax- 
ies at fixed a are the younger, Fe-rich galaxies while the 
fainter galaxies are older and more Fe-poor. This con- 
firms that the observed anti-correlation of age and [Fe/H] 
must be genuine, rather than an effect of correlated stel- 
lar population modelling errors. 



The CO- variation of age, [Fe/H], and [IVIg/Fe] we ob- 
serve at fixed cr cannot be driven by large-scale environ- 
ment. The age-[Fe/H] anti-correlation is observed sepa- 
rately in the Trager et al. (2000) data, which include only 
field galaxies, and in the Smith et al. (2007b) data, which 
include only cluster galaxies. It is also observed in our 
stacked spectra, which average over all environments. 

The observed co- variation of age, [Fe/H], and []VIg/Fe] 
is qualitatively consistent with a scenario in which all 
galaxies at the same a start forming stars at the same 
z but the duration of star formation varies between 
galaxies. At fixed cr, galaxies with short star forma- 
tion timescales end star formation early and consequently 
have older mean stellar ages. They also have less time for 
enrichment in SN la products, resulting in lower [Fe/H] 
and higher [JVIg/Fe]. To the extent that a measures the 
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depth of the gravitational potential, the loss of metals 
through winds should be consistent for all galaxies at 
the same a if they experience comparable galactic winds. 
Thus early truncation of star formation should have little 
effect on [Mg/H] but should result in a genuine under- 
abundance of [Fe/H], as observed. To also be consis- 
tent with the general trends of Figure 6, these variations 
in the duration of star formation must be superimposed 
on a staged star formation model, similar to the Noeske 
et al. (2007) model, where the redshift at which star for- 
mation peaks and the duration of star formation depend 
on a. Unlike the Noeske model, the data presented here 
additionally require a spread in star formation duration 
at fixed (7. A quantitative test of such a model requires 
the construction of chemical evolution models, which we 
defer to future work. Furthermore, this scenario postu- 
lates a range of star formation timescales at fixed a but 
does not propose a physical mechanism to produce the 
observed variation. 

A key result of this analysis is that the co- variation of 
age, [Fe/H], and [Mg/Fe] at fixed a also correlates with 
galaxy luminosity, such that the brighter galaxies are 
younger, Fe-rich, and less a-enhanced than the fainter 
galaxies at the same u. If the above scenario is correct, 
the more luminous galaxies at fixed a would be those 
which experienced more extended star formation histo- 
ries. This at first appears to be a trivial result, since 
galaxies with more extended star formation and younger 
mean ages are expected to have lower MjL and therefore 
higher L at fixed a. However, the differences in MjL 
predicted from stellar population models are too small 
to account for the observed variation in L. This will be 
shown in detail in future papers in this series, but can 
also be inferred from Figure 3. This figure shows that 
the range in galaxy colors at fixed u is small, suggest- 
ing that the range of MjL is likewise small. In fact, 
the variation in L at fixed a turns out to be dominated 
by variation in total stellar mass at fixed cr, while MjL 
variations contribute only a modest amount. Thus any 
physical mechanism invoked to explain the observed co- 
variation in stellar population properties at fixed a must 
also reproduce differing total stellar mass content in the 
galaxies. 

7. CONCLUSIONS 

This analysis has used very high SjN stacked spec- 
tra of ~ 16,000 SDSS early type galaxies to map out 
variations in stellar population age, [Fe/H], [Mg/H], and 
[Mg/Fe] with galaxy cr, L, and color. In addition to 
stiidying the mean trends in stellar populations with 
each of these! three global properties individually, we 
have mapped out variations in age, [Fe/H], [Mg/H], and 
[Mg/Fe] in the three dimensional cr-L-color parameter 
space. This allows us to understand the differing behav- 
ior of stellar populations as a function of a versus L, and 
to explore the systematic variations in stellar populations 
at fixed a. 

We find the following results for quiescent galaxies: 

1. Luminosity and a both increase with galaxy mass, 
yet these "size" measures are not the same. Fun- 
damental stellar population variables such as age, 
[Fe/H], [Mg/H], and [Mg/Fe] scale differently with 
respect to L and to a, and also with respect to 
galaxy color. 



2. Higher cr galaxies are typically more luminous and 
redder than lower cr galaxies. However, at fixed cr, 
there is no color-magnitude relation, despite the 
substantial variation in M^. The color-magnitude 
relation of passive galaxies is therefore a result of 
combining the a-L and cr-color relations, in agree- 
ment with Bernardi et al. (2005). 

3. Mean luminosity-weighted stellar population age, 
[Fe/H], [Mg/H], and [Mg/Fe] aU correlate with a 
such that galaxies with higher cr tend to be older 
and have higher [Fe/H], [Mg/H], and [Mg/Fe]. At 
fixed cr, the spread in [Mg/H] and [Mg/Fe] is small, 
while both [Fe/H] and age show significant spread. 

4. At fixed cr, brighter galaxies have lower age, 
higher [Fe/H], slightly higher [Mg/H], and lower 
[Mg/Fe] than their fainter counterparts. The anti- 
correlation of age and [Fe/H] conspires to keep the 
color-magnitude relation flat at fixed u and con- 
tributes to the overall narrowness of the color-a 
relation. 

5. At fixed cr, age is also correlated with color such 
that younger galaxies are both brighter and bluer 
at fixed cr than the fainter, redder, older galaxies. 
In contrast, [Fe/H] and [Mg/H] do not vary sys- 
tematically with color, while [Mg/Fe] varies only 
mildly with color. 

6. The variation in stellar population properties with 
L at fixed cr acts to reinforce the cr [Fe/H] trend, 
resulting in a strong and tight L-[Fe/H] correla- 
tion. In contrast, the u-age and cr- [Mg/Fe] cor- 
relations are opposed by the variations with L at 
fixed cr, leaving only weak L-age and L- [Mg/Fe] 
correlations. In a similar way, residuals in color cor- 
relate with stellar population properties such that 
there is a strong color-age trend but only weak 
color-[Fe/H] and color-[Mg/Fe] trends. Age corre- 
lates most closely with color, [Fe/H] correlates most 
closely with L, and [Mg/H] and [Mg/Fe] correlate 
most closely with a. Only u (and not L or color) 
correlates strongly with all four stellar population 
properties. 

7. Trends in age, [Fe/H], [Mg/H], and [Mg/Fe] at fixed 
cr are likely not driven by environment. The age- 
[Fe/H] anti-correlation at fixed cr is detected in 
samples of individual galaxies which reside in sim- 
ilar environments (e.g., the field galaxies in Trager 
et al. 2000 and the cluster galaxies in Smith et al. 
2007b) and also in our sample of stacked spectra, 
which average over all environments. 

The variations in stellar population properties at fixed 
cr presented here illustrate that the narrow, seemingly 
one-dimensional color magnitude relation of quiescent 

galaxies conceals an underlying set of star formation his- 
tories that populate a two-parameter family. Not only 
do the star formation histories of galaxies vary system- 
atically with their cr, but there is clearly a range of pro- 
cesses at work in galaxies of the same cr. These result 
in an anti-correlation between galaxy age and [Fe/H] at 
fixed o and a weaker but positive age- [Mg/Fe] correla- 
tion. Moreover, these properties are linked to the present 
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day luminosities of the galaxies, such that brighter galax- 
ies arc younger, more metal-rich, and less a-enhanced 
than fainter galaxies at the same a. The companion pa- 
pers in this series will explore in more detail the con- 
nection between galaxy structure and galaxy stellar pop- 
ulations, providing evidence for the ways in which the 
star formation history of a galaxy is linked to its mass 
assembly and structural evolution. 

The authors would like to thank Rcnbin Yan for pro- 
viding the emission line measurements used to identify 
the sample of early-type galaxies used here, as well as 
an anonymous referee for provided valuable feedback. 
This work was supported by National Science Founda- 
tion grant AST 05-07483. 

Funding for the creation and distribution of the SDSS 
Archive has been provided by the Aired P. Sloan Foun- 



Adelman-McCarthy, J. K., ct al. 2006, ApJS, 162, 38 
Adelman-McCarthy, J. K., ct al. 2007, ArXiv c-prints, 707 
Baldwin, J. A., Phillips, M. M., & Tcrlcvich, R. 1981, PASP, 93, 5 
Bender, R., Burstcin, D., & Fabcr, S. M. 1993, ApJ, 411, 153 
Bernardi, M. 2007, AJ, 133, 1954 

Bernardi, M., Hyde, J. B., Shcth, R. K., Miller, C. J., & Nichol, 

R. C. 2007, AJ, 133, 1741 
Bernardi, M., et al. 2003a, AJ, 125, 1817 
Bernardi, M., et al. 2003b, AJ, 125, 1882 

Bernardi, M., Slieth, R. K., Nichol, R. C, Schneider, D. P., & 

Brinkmann, J. 2005, AJ, 129, 61 
Blanton, M. R., ct al. 2003, AJ, 125, 2348 
Blanton, M. R., ct al. 2001, AJ, 121. 2358 
Blanton, M. R., ct al. 2005, AJ, 129, 2562 

Bower, R. G., Luccy, J. R., & EUis, R. S. 1992, MNRAS, 254, 601 
Burbidge, E. M., Burbidge, G. R., & Fish, R. A. 1961, ApJ, 133, 
393 

Burstein, D., Faber, S. M., Gaskell, C. M., & Krumm, N. 1984, 
ApJ, 287, 586 

Cardial, N., Gorgas, J., Cenarro, J., & Gonzalez, J. J. 1998, A&AS, 
127, 597 

CoUess, M., Burstein, D., Davies, R. L., McMahan, R. K., Saglia, 

R. P., & Wegner, G. 1999, MNRAS, 303, 813 
Cool, R. J., Eisenstein, D. J., Johnston, D., Scranton, R., 

Brinkmann, J., Schneider, D. P., k. Zehavi, I. 2006, AJ, 131, 

736 

Djorgovski, S. & Davis, M. 1987, ApJ, 313, 59 

Dressier, A., Lyndcn-BcU, D., Burstcin, D., Davies, R. L., Faber, 

S. M., Terlevich, R., & Wcgncr, G. 1987, ApJ, 313, 42 
Faber, S. M. 1973, ApJ, 179, 731 

Faber, S. M., Dressier, A., Davies, R. L., Burstein, D., & Lynden- 

Bell, D. 1987, in Nearly Normal Galaxies; From the Planck Time 

to the Present, ed. S. M. Fabcr, 175-183 
Fabcr, S. M. & Jackson, R. E. 1976, ApJ, 204, 668 
Fabcr, S. M., Trager, S. C., Gonzalez, J. J., & Worthey, G. 1995, 

in lAU Symposium, Vol. 164, Stellar Populations, ed. P. C. van 

dcr Kruit & G. Gilmorc, 249 
Fukugita, M., Icliikawa, T., Gunn, J. E., Doi, M., ShimasaJfu, K., 

& Schneider, D. P. 1996, AJ, 111, 1748 
Pulbright, J. P., McWiUiam, A., & Rich, R. M. 2007, ApJ, 661, 

1152 

Gallazzi, A., Chariot, S., Brinchmann, J., & White, S. D. M. 2006, 

MNRAS, 370, 1106 
Gallazzi, A., Chariot, S., Brinchmann, J., White, S. D. M., & 

Tremonti, C. A. 2005, MNRAS, 362, 41 
Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&AS, 

141, 371 

Gladders, M. D. & Yce, H. K. C. 2000, AJ, 120, 2148 
Graves, G. J., Faber, S. M., Schiavon, R. P., & Yan, R. 2007, ApJ, 
671, 243 

Graves, G. J. & Schiavon, R. P. 2008, ArXiv e-prints, 802 
Gunn, J. E., et al. 1998, AJ, 116, 3040 
Gunn, J. E., et al. 2006, AJ, 131, 2332 



dation, the Participating Institutions, the National Aero- 
nautics and Space Administration, the National Science 
Foundation, the US Department of Energy, the Japanese 
Monbukagakusho, and the Max-Planck Society. The 
SDSS Web site is http://www.sdss.org/. 

The SDSS is managed by the Astrophysical Research 
Consortium (ARC) for the Participating Institutions. 
The Participating Institutions are the University of 
Chicago, Fermilab, the Institute for Advanced Study, 
the Japan Participation Group, the Johns Hopkins Uni- 
versity, the Korean Scientist Group, Los Alamos Na- 
tional Laboratory, the Max-Planck-Institute for Astron- 
omy (MPIA) , the Max-Planck-Institute for Astrophysics 
(MPA), New Mexico State University, University of 
Pittsburgh, University of Portsmouth, Princeton Uni- 
versity, the United States Naval Observatory, and the 
University of Washington. 



Hogg, D. W., Finkbciner, D. P., Schlegel, D. J., & Gunn, J. E. 

2001, AJ, 122, 2129 
Humphrey, P. J. & Buotc, D. A. 2006, ApJ, 639, 136 
Ivczic, Z., ct al. 2004, Astronomische Nachrichten, 325, 583 
Jones, C. D. 1999, PhD thesis, Univ. of Washington 
J0rgensen, I. 1999, MNRAS, 306, 607 

J0rgensen, I., Pranx, M., & Kjaergaard, P. 1995, MNRAS, 276, 
1341 

Kauffmann, G., ct al. 2003, MNRAS, 346, 1055 

Kelson, D. D., lUingworth, G. D., Franx, M., & van Dokkum, P. G. 

2006, ApJ, 653, 159 
Kodama, T. & Arimoto, N. 1997, A&A, 320, 41 
Kormcndy, J. 1985, ApJ, 295, 73 

Korn, A. J., Maraston, C, & Thomas, D. 2005, A&A, 438, 685 
Kuntschncr, H. & Davies, R. L. 1998, MNRAS, 295, L29 
Kuntschner, H., Lucey, J. R., Smith, R. J., Hudson, M. J., & 

Davies, R. L. 2001, MNRAS, 323, 615 
Lauer, T. R. 1985, ApJ, 292, 104 
Lauer, T. R., ct al. 2007, ApJ, 662, 808 

Lisker, T., Grcbcl, E. K., BinggeU, B., & Glatt, K. 2007, ApJ, 660, 
1186 

Lupton, R., Gunn, J. E., Ivezic, Z., Knapp, G. R., & Kent, S. 

2001, in Astro. Soc. of the Pacific Conf. Ser., Vol. 238, ed. F. R. 

Harnden, Jr., F. A. Primini, & H. E. Payne, 269 
Nelan, J. E., Smith, R. J., Hudson, M. J., Wegner, G. A., Lucey, 

J. R., Moore, S. A. W., Quinney, S. J., & Suntzeff, N. B. 2005, 

ApJ, 6,32. 137 
Nocskc, K. G., ct al. 2007, ApJ, 660, L47 

Pier, J. R., Munn, J. A., Hindslcy, R. B., Hennessy, G. S., Kent, 

S. M., Lupton, R. H., & Ivezic, Z. 2003, AJ, 125, 1559 
Poggianti, B. M., et al. 2001, ApJ, 562, 689 
Rix, H.-W. & White, S. D. M. 1992, MNRAS, 254, 389 
Sandage, A. & Visvanathan, N. 1978, ApJ, 223, 707 
Schawinski, K., Thomas, D., Sarzi, M., Maraston, C, Kaviraj, S., 

Joo, S.-J., Yi, S. K., & Silk, J. 2007, MNRAS, 382, 1415 
Schiavon, R. P. 2007, ApJS, 171, 146 
Smith, J. A., et al. 2002, AJ, 123, 2121 

Smith, R. J., Lucey, J. R., & Hudson, M. J. 2007a, MNRAS, 381, 
1035 

— . 2007b, ArXiv e-prints, 712 

Stoughton, C, et al. 2002, Proceedings of the SPIE Vol. 4836, ed. 

Tyson, J. A.; & Wolff, S., 339-349 
Strauss, M. A., et al. 2002, AJ, 124, 1810 
Tcrlcvich, A. I. & Forbes, D. A. 2002, MNRAS, 330, 547 
Thomas, D., Maraston, C, Bender, R., & Mendes de Oliveira, C. 

2005, ApJ, 621, 673 
Tinsley, B. M. 1981, MNRAS, 194, 63 

Trager, S. C, Faber, S. M., & Dressier, A. 2008, MNRAS, 441 
Trager, S. C, Faber, S. M., Worthey, G., & Gonzalez, J. J. 2000, 
AJ, 120, 165 

Trager, S. C, Worthey, G., Faber, S. M., Burstein, D., & Gonzalez, 

J. J. 1998, ApJS, 116, 1 
Tucker, D. L., et al. 2006, Astronomische Nachrichten, 327, 821 
Worthey, G. 1994, ApJS, 95, 107 



20 



Graves et al. 



8000 
6000 

4000 h 

eooo - 



fjL - 0.001 



u = 0,0069 




400 



-i 100 





-0,04 -0,02 0.00 0,0S 0,04 -0.02 0,00 0.02 0,04 



Fig. Al. — Distribution of g — r color differences between composite bulge + disk and de Vaucouleurs model light profiles for our sample 
galaxies. In the bottom panel, the galaxies with color differences < 0.001 mag have been removed, showing the nearly Gaussian distribution 
of the remaining color differences. Overall, > 97% of the color differences are less than the 0.022 mag rms error in the g — r colors measured 
from SDSS photometry (dotted lines). Wc conclude that using magnitudes and colors derived from de Vaucouleurs fits to the galaxy light 
profiles does not substantially affect the results presented here. 
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APPENDIX 



BULGE + DISK VS. DE VAUCOULEURS MAGNITUDES AND COLORS 

This analysis is based on magnitudes and colors measured from de Vaucouleur fits to the galaxy light profiles. 
However, given that a fraction of the galaxies in the sample have a visible disk component it is important to verify 
that the de Vaucouleurs magnitudes and colors are a reasonable approximation to the total galaxy magnitudes and 
colors. 

The Petrosian magnitudes measured by the SDSS pipeline have larger errors (typically 5% in the r-band) than the de 
Vaucouleurs model fits (typically 2% in the r-band) as determined by repeat observations of targeted galaxies. However, 
they have the benefit of being independent of the galaxy light profile. We compare Petrosian r-band magnitudes to 
the de Vaucouleurs magnitudes used in our analysis and find that the dc Vaucouleurs magnitudes are systematically 
brighter by < 0.1 mag (consistent with Blanton et al. 2001), with rms scatter of 0.04 mag. This scatter is much 
smaller than the 0.8 mag bin widths used to sort and stack galaxy spectra and should therefore have little impact on 
the analysis presented here. 

Color effects may be more important, both because the color differences observed in our sample galaxies are intrinsi- 
cally small and because disk galaxies may exhibit substantial color gradients. Wc combine the separate dc Vaucouleurs 
and exponential fits to galaxy light profiles, weighting each component by the fractional contribution of each com- 
ponent (the SDSS pipeline value fracDeV) to produce composite bulge -|- disk galaxy magnitudes and colors. We 
compare the g — r colors from the bulge -|- disk models to the de Vaucouleurs profile colors in Figure Al. In the 
top panel we show the distribution of color differences for all galaxies that go into our stacked spectra (i.e., after 
the color- and magnitude-outliers have been removed from each a bin). The majority of galaxies show < 0.001 mag 
difference between the composite and dc Vaucouleur color values because the galaxies are highly bulge-dominated. 
Removing these and examining the distribution of color difFcrcnccs for galaxies with color differences larger than 0.001 
mag (bottom panel), we find that the distribution is roughly gaussian, with ji = 0.0015 mag and cr = 0.0069 mag. 

Given that the rms error for 5 — r color photometry in SDSS is ~ 2% (0.022 mag, dotted lines in Figure Al), the 
color differences from the different model light profiles are within this error for > 97% of our galaxy sample. We 
therefore conclude that the use of de Vaucouleurs model fits for the magnitudes and colors in this analysis does not 
have a significant effect on the results. 
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Fig. B1. — Stellar population model grid showing the effect of adding a young burst to an underlying old population. The model grid 
is for a solar-scale abundance pattern model from Schiavon (2007). Solid (dotted) lines connect models with the same [Fe/H] (age). The 
gray track shows the effect of taking an old population with age = 14.1 Gyr and [Fe/H] = -0.4, and adding to it varying fractions of a 
1.2 Gyr burst with the same [Fe/H]. Gray diamonds indicate the fraction of the total luminosity at 5000A provided by the young burst, 
starting at 0% at the bottom and increasing in increments of 5%. Index values from the stacked spectra are overplotted as triangles, with 
colors indicating the u range of the stacked spectrum, as in Figure 3 (purple shows the lowest a bin, red the highest cr bin). The three 
data points at fixed u show the three different bins in L. The age-[Fe/H] anticorrelation observed at fixed a cannot simply be due to the 
presence of a young sub-population in the brighter galaxies because the trajectory of a young burst at constant [Fe/H] is nearly orthogonal 
to the observed co-variation of (Fe) and H/3 at fixed cr. The brighter galaxies must have higher average [Fe/H] than their faint counterparts 
at the same cr, in addition to having younger mean ages. 



THE EFFECT OF YOUNG SUB-POPULATIONS 

Because the observed age-[Fe/H] anti-correlation at fixed cr is in the same sense as the age-metalhcity degeneracies 
that plague stellar population studies, it is necessary to confirm that systematic errors in the modelling process are 
not producing a spurious anti-correlation. We showed in §13 that the age-[Fe/H] anti-correlation at fixed a is not due 
to random observational errors. However, the models used in this work are single stellar population (SSP) models, 
whereas galaxies are composite populations with extended star formation histories. The mean age measured from 
an SSP model is a mean luminosity- weighted age for all the stars contributing to the total light. Because young 
populations have lower mass-to-light ratios and stronger Balmer absorption lines, a sub-population of intermediate 
age stars (either a younger sub-population within individual galaxies or a subset of significantly younger galaxies in 
the bin) will contribute disproportionately to the total galaxy light, skewing the SSP age to low values. If the young 
sub-population greatly strengthens the Balmer absorption lines without changing the Fe absorption lines much, the 
resulting model grid inversion may infer higher metallicities just due to the stronger Balmer lines rather than a genuine 
difference in metallicity. 

This is in the correct sense to account for the anti-correlation we see between galaxy age and [Fe/H] at fixed a. 
However, it alone cannot account for metallicity differences at fixed a, as illustrated in Figure Bl. This figure shows a 
stellar population model grid for the (Fe) and H/3 indices. Solid lines connect models of constant [Fc/H], while dotted 
lines connect models of constant age. Triangles show the index measurements for stacked spectra in this work, with 
colors indicating the a range represented (purple shows the lowest a bin, red the highest a bin, as in Figure 3). We 
have averaged over the three different color bins so the three points shown for each a bin represent three different 
sub-bins in L. The gray track with diamonds shows the effect of adding a young sub-population to an underlying old 
population. The old population has age — 14.1 Gyr and [Fe/H] — -0.4, while the young sub-population has the same 
[Fe/H] and age = 1.2 Gyr. Gray diamonds indicate the fraction of the total luminosity at 5000 A provided by the 
burst, starting at 0% at the bottom and increasing in increments of 5%. 

It is clear that the trajectory in index-index space covered by composite burst models with fixed [Fe/H] is very 
different from the observed co- variation in galaxy age and [Fe/H]. It is true that the [Fe/H] inferred for a composite 
burst would be skewed slightly high, since the gray track in Figure Bl lies consistently to the right of the [Fe/H] = 
-0.4 grid line despite the fact that both sub-populations in the composite models have [Fe/H] = -0.4. However, the 
effect is small compared to the observed variation in (Fe) at fixed a in the galaxy sample. The data therefore require 
genuine variations in mean [Fe/H] at fixed cr, in addition to age variations. 
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g — r color as measured from the total dc Vaucouleurs galaxy photometry 
Total number of galaxies in the staekcd spectrum. 
^Effective median S/N of the stacked spectrum. 
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TABLE C2 

Lick Indices and Stellar Population Properties of Stacked Spectra 



H/3 


<Fe) 


Mg b 


Age 


[Fe/H] 


[Mg/H] 


[Mg/Fe] 


(A) 


(A) 


(A) 


(Gyr) 


(dex) 


(dex) 


(dex) 





2 


1 n 

. J-U 


-1- 


u.uo 


2 


.o^ 


-1- 


OR 

u.uo 


0. lU 


-1- 


Of> 
u.uo 


0.0 


-1- 


^ 
u.o 


-0 22 


-1- 


0*^ 

u.uo 


n 
-u. 


.uo 




OA 
u.u^ 


n 

U. 


14 


-i- 


0^ 
U.UO 


2 


2 


.09 


lb 


0.04 


2, 


.42 


-|- 


0.05 


3.38 


-|- 


0.04 


5.5 


d= 


0.4 


-0.16 


J- 


0.03 


-0. 


.01 


dz 


0.03 


0. 


.15 


-j- 


0.02 


3 


]^ 


.88 


-|- 


0.06 


2, 


.58 


-|- 


0.07 


3.78 


-|- 


0.06 


8.0 


-|- 


0.8 


-0.12 


-|- 


0.03 


0. 


.03 


-|- 


0.05 


0. 


15 


-|- 


0.03 




2 


1 7 


-1- 


0*^ 


2 


.oy 




OA 


9 QQ 
z.yy 


-1- 


0*^ 
u.uo 


^ 
o.u 




'i 
u.o 


1 
-u. 10 




09 
u.uz 


n 

-u. 




-i_ 

3Z 


0'^ 

u.uo 


n 
u. 


07 
U ( 


j- 

3Z 


09 
u.uz 


c; 
o 


2 


07 




0*^ 

u.uo 


2 


.oo 




0*^ 

u.uo 


'i 90 
o.zu 


-1- 


0*^ 

u.uo 


R 

o.u 


-1- 


A 


-0 24 




n 09 
u.uz 


n 

-u. 


1 n 


J- 


O'^ 

u.uo 


n 

u. 


14 


-L 

zn 


09 
u.uz 


^; 
u 


^ 


.yu 


-1- 


0*^ 

u.uo 


2 


.'iO 




OA 


'i Al 


-1- 


0*^ 

u.uo 


7 1 
1 . 1 


-1- 


A 


-0 21 


31 


n 09 

u.uz 


n 

-u. 


.U 1 


-1- 


O'? 

u.uo 


n 

u. 


14 


-1- 


09 
U.uz 






.98 




0.06 


2 


.24 


=b 


0.07 


2.88 


-|- 


0.06 


7.5 


-|- 


0.9 


-0.34 




04 


-0 


.25 


zb 


0.05 


0. 


09 


zb 


0.04 






89 




04 


2 


23 




0.05 


3.14 


zt 


0.04 


8.7 




7 


-0 38 




02 


-0 


.22 


zb 


0.03 


0. 


16 


zb 


0.02 


9 




so 




nfi 
u.uu 


2 


•^0 

. ou 


-1- 


n OR 
u.uu 


A1 




0^^ 
u.uo 


lU.O 


-j- 


1 1 


-u.oo 


-1- 


n 0*^ 
u.uo 


n 

-u. 


1 7 




0^^ 

u.uo 


n 

u. 


1 s 

,10 


-1- 


0^^ 
u.uo 


10 


2 


.16 


± 


0.04 


2, 


.49 


-|- 


0.04 


3.34 


-|- 


0.04 


4.8 


-|- 


0.3 


-0.10 


-|- 


0.02 


0. 


02 




0.03 


0. 


,12 


-|- 


0.02 


11 


1 


.92 


-j- 


0.03 


2, 


.51 


-|- 


0.03 


3.48 


-|- 


0.03 


7.5 


-|- 


0.4 


-0.17 


=b 


0.02 


-0 


.04 




0.03 


0. 


13 


zb 


0.02 


12 


\ 


.77 


-l- 


0.04 


2, 


.70 


-|- 


0.04 


3.92 


-|- 


0.04 


9.6 




0.7 


-0.09 


=b 


0.04 


0. 


03 




0.05 


0. 


12 




0.02 




2 


07 




u.uo 


2 


. oo 




0*^ 
u.uo 


0.10 


-1- 


0*^ 
u.uo 


o.y 




A 


-0 24 




09 
u.uz 


n 

-u. 


1 n 




0'^ 
u.uo 


n 
u. 


14 




09 
U.UZ 


14 


]^ 


.yo 


-1- 


09 


2 


47 


-1- 


0*^ 

u.uo 


3 44 


-1- 
IC 


0*^ 

u.uo 


7 1 


-1- 




u.o 


01s 
-u. 10 


-1- 


01 
U.Ul 


n 
-u. 


.uo 


J- 


09 
U.UZ 


n 

U. 


1 'i 
. 10 


J- 
^ 


01 

U.Ul 


15 


]^ 


.82 


-l- 


0.03 


2. 


.47 


-|- 


0.03 


3.61 


-|- 


0.03 


9.1 


-|- 


0.5 


-0.22 


d= 


0.02 


-0. 


.05 


dz 


0.03 


0. 


.17 


-j- 


0.02 


16 


1 


.85 


-l- 


0.05 


2, 


.33 


-|- 


0.05 


3.10 


-|- 


0.04 


9.3 


-|- 


0.8 


-0.33 


-U 

IE 


0.03 


-0. 


.24 


zb 


0.04 


0. 


.09 


-|- 


0.03 


17 


\ 


.85 


-l- 


0.03 


2, 


.30 


-|- 


0.04 


3.42 


-|- 


0.03 


9.3 


-|- 


0.6 


-0.33 


-|- 


0.02 


-0. 


.13 


-j- 


0.03 


0. 


.20 


-j- 


0.02 


1 s 




77 




u.uo 


2 


'^7 

.O 1 


4- 


u.uo 


o.ou 


-1- 


Of> 
u.uo 


10 ^ 


4- 


Q 
u.y 


'^1 

-U.Ol 


4- 


0^ 
u.uo 


n 
-u. 


]^]^ 




OA 


U. 


90 

.ZU 


-L. 
^ 


0*^ 
u.uo 


1 Q 


2 


.uo 


-L. 


09 
u.uz 


2 


.oo 


-1- 


09 
u.uz 


^ RO 
o.ou 


-1- 


09 
u.uz 


^ R 
0.0 


-1- 
-L 




u.o 


OS 

-U.UO 


-1- 


09 
u.uz 


u. 


C\f{ 
uo 


-1- 


09 
u.uz 




u. 


14 


-1- 


01 

U.Ul 


90 




.oo 


-1- 


09 
u.uz 


2 




-1- 


09 
u.uz 


o.ou 


-1- 


09 
u.uz 


S A 


-L. 




u.o 


17 
-u. 1 / 




01 
U.Ul 


-0, 


.05 


-1- 

3Z 


09 
U.uz 


A 
U. 


12 


-1- 


01 
U.Ul 


21 


^ 


7'^ 


-i~ 
3Z 


0*^ 

u.uo 


2 


70 

. / u 




0*^ 

u.uo 


'i Q7 
o.y ( 




0*^ 

u.uo 


10 1 
lU. 1 




^ 

u.o 


Oil 

-U. 1 1 




01 
U.Ul 


0, 


,03 


=h 


0*^ 

U.UO 


n 

u. 


14 


-1- 
3Z 


09 
U.UZ 


22 


2 


on 
.uu 


J- 


09 
u.uz 


2 


AQ 




09 
u.uz 


o.oy 




09 
u.uz 


R R 

u.u 






u.o 


-0 14 


31 


n 01 

U.Ul 


0, 


,02 


zc 


09 
U.UZ 


n 
u. 


1 R 

. lU 


4- 

31 


09 
U.UZ 


23 


]^ 


.86 


-|- 


0.02 


2, 


.50 


-|- 


0.02 


3.67 


-|- 


0.02 


8.4 


-|- 


0.3 


-0.18 


-|- 


0.01 


-0, 


.01 


± 


0.01 


0. 


17 


dz 


0.01 


24 


]^ 


.75 


-|- 


0.02 


2, 


.65 


-|- 


0.02 


3.92 


-|- 


0.02 


9.9 


-|- 


0.4 


-0.13 


-|- 


0.01 


0, 


,02 


-|- 


0.02 


0. 


,15 


-|- 


0.02 




]^ 


S7 




0"^ 

u.uo 


2 


9Q 
.zy 




OA 


o.oy 




0"^ 
u.uo 


s 7 




R 
u.u 


'^9 
-u.oz 




09 
u.uz 


-0, 


.12 


-1- 


0*^ 
u.uo 


n 

u. 


90 
,zu 


-1- 


09 
u.uz 


9fi 




S9 




0"^ 

u.uo 


2 


f^7 
.0 i 




0"^ 

u.uo 


0.00 




0"^ 

u.uo 


Q 

y .u 




A 

u.^ 


1 
-u. 10 




01 
U.Ul 


0, 


,01 


-1- 


0*^ 

u.uo 


u. 


1 R 

. lU 


zn 


09 
U.uz 


97 




.UO 




0"^ 

u.uo 


2 


49 




OA 


Q7 
O.y ( 




0"^ 

u.uo 


12 3 




7 
u. / 


9R 
-u.zu 




n 09 

u.uz 


-0, 


.03 


j- 

ZjZ 


0*^ 

u.uo 


u. 


9*^ 

. zo 


-1- 

3Z 


09 
U.uz 


28 




.95 


-|- 


0.02 


2, 


.61 


lb 


0.02 


3.84 


-|- 


0.02 


6.9 


-|- 


0.2 


-0.06 


-|- 


0.01 


0, 


,11 


-|- 


0.02 


0. 


17 


-|- 


0.02 


9Q 


]^ 


7Q 

. i y 




01 
U.Ul 


2 


.uo 




09 

u.uz 


SA 




01 
U.Ul 


Q 9 
y . z 




9 
u.z 


-0 12 


31 


n 01 

U.Ul 


0, 


,02 


-1- 


09 
u.uz 


u. 


14 


-1- 


09 
u.uz 


ou 


]^ 


79 




09 
u.uz 


2 


.uo 




09 

u.uz 


4 10 




09 
u.uz 


1 

lU.O 




A 


oil 
-u. 1 1 




01 
U.Ul 


0, 


.07 


-1- 


09 
U.UZ 


u. 


1 s 

.10 


-1- 


09 
U.uz 


Ol 


]^ 


Q1 

.yi 


-1- 


09 

u.uz 


2 


.00 


-1- 


09 

u.uz 


'i 7R 
0. ( 


-1- 


09 

u.uz 


7 ^ 
1 .0 


-1- 


9 
U.z 


-0 14 


-1- 


01 
U.Ul 


0, 


.04 




09 
U.uz 


n 

u. 


1 s 

.10 




09 
U.uz 




]^ 


7Q 
. / y 


-1- 


01 
U.Ul 


2 


.ou 


-1- 


01 
U.Ul 


A 0"^ 


-1- 


01 
U.Ul 


Q 1 

y . 1 


-1- 


9 
U.z 


1 
-u. 10 


-1- 


01 
U.Ul 


0, 


.08 




01 
U.Ul 


n 
u. 


21 




01 
U.Ul 




]^ 


.70 


-|- 


0.01 


2, 


.69 


-|- 


0.02 


4.23 


-|- 


0.01 


10.6 


-|- 


0.3 


-0.10 


-|- 


0.01 


0. 


.10 


-j- 


0.02 


0. 


.20 


-j- 


0.02 




]^ 


78 
. / o 


-1- 


0^ 
u.uo 


2 


.00 


-1- 


0^ 
u.uo 


^ 89 


-1- 


0^ 
u.uo 


10 

lU.U 


4- 


^ 
u.o 


98 
-u.zo 


J- 


01 

U.Ul 


-0.03 


^ 


09 
u.uz 


n 

u. 


9^ 

.zo 


^ 


09 
u.uz 


35 


]^ 


.68 


lb 


0.02 


2, 


.49 


-|- 


0.02 


4.05 


-|- 


0.02 


11.5 


J- 


0.4 


-0.24 


-|- 


0.01 


0. 


.02 


-|- 


0.02 


0. 


.26 


-|- 


0.02 


ou 




R7 


-L 
31 


0^ 
u.uo 


2 


fiO 

.ou 


-1- 


0^ 
u.uo 


4 28 


-1- 


0^ 
u.uo 


11^ 
11.0 


-L 
31 


^ 
u.o 


17 
-u.i 1 


-1- 


09 
U.UZ 


0. 


.10 


-1- 


0"^ 
u.uo 


u. 


97 
.z / 


-1- 


09 
u.uz 


o / 




8R 
.Oo 




09 


2 


.00 


-1- 


09 
u.uz 


4 16 


-1- 


09 
u.uz 


7 7 


4- 


9 
U.z 


0^ 

-u.uo 


-1- 


01 

U.Ul 


0. 


.17 


-I- 


01 

U.Ul 


n 
u. 


90 
.ZU 


-L. 


01 

U.Ul 


oo 




71 




01 

U.Ul 


2 


71 

. / 1 


-1- 


09 
u.uz 


A 97 


-1- 


01 

U.Ul 


10 2 


4- 


9 
U.z 


0>^ 

-u.uo 


-1- 


01 

U.Ul 


0. 


.13 


-i- 


01 

U.Ul 


u. 


21 




01 

U.Ul 


oy 




RO 


-1- 


09 
u.uz 


2 


71 

. i 1 


-1- 


09 
u.uz 


4 42 


-1- 


09 
u.uz 


12 3 


-1- 
^ 


^ 

u.o 


-0 12 


-1- 


01 
U.Ul 


0, 


.12 


-1- 


09 
U.UZ 


n 

u. 


24 


-1- 


01 
U.Ul 


40 


1 


.80 




0.02 


2 


.65 


-|- 


0.02 


4.14 


-|- 


0.02 


8.8 


J- 


0.2 


-0.08 


-|- 


0.01 


0, 


.13 


-|- 


0.01 


0. 


.21 


-|- 


0.01 


41 


1, 


.69 


ib 


0.01 


2, 


.69 


ib 


0.01 


4.31 


ib 


0.01 


10.6 


ib 


0.2 


-0.10 


zb 


0.01 


0, 


.12 


zb 


0.02 


0, 


.22 


zb 


0.02 


42 


1, 


.62 


lb 


0.02 


2, 


.72 


lb 


0.02 


4.46 


lb 


0.02 


11.8 


lb 


0.3 


-0.10 


zb 


0.01 


0, 


.14 


zb 


0.01 


0, 


.24 


zb 


0.01 


43 


1, 


.69 


lb 


0.03 


2, 


.45 


lb 


0.03 


4.19 


lb 


0.03 


11.4 


lb 


0.6 


-0.25 


zb 


0.02 


0, 


.06 


zb 


0.03 


0, 


.31 


zb 


0.02 


44 


1, 


.64 


lb 


0.02 


2, 


.58 


lb 


0.02 


4.33 


lb 


0.02 


11.9 


lb 


0.4 


-0.18 


zb 


0.01 


0, 


.10 


zb 


0.01 


0, 


.28 


zb 


0.01 


45 


1, 


.57 


=b 


0.02 


2, 


.67 


lb 


0.03 


4.54 


lb 


0.02 


13.1 


lb 


0.6 


-0.14 


zb 


0.01 


0, 


.13 


zb 


0.02 


0, 


.27 


zb 


0.02 


46 


1, 


.67 


lb 


0.03 


2, 


.73 


lb 


0.03 


4.48 


lb 


0.03 


10.9 


lb 


0.6 


-0.08 


zb 


0.02 


0, 


.17 


zb 


0.03 


0, 


.25 


zb 


0.02 


47 


1, 


.65 


± 


0.02 


2, 


.78 


lb 


0.02 


4.48 


lb 


0.02 


10.9 


zb 


0.4 


-0.05 


zb 


0.01 


0, 


.19 


zb 


0.02 


0, 


.23 


zb 


0.02 


48 


1, 


.59 


± 


0.03 


2, 


.74 


lb 


0.03 


4.69 


lb 


0.03 


12.2 


lb 


0.6 


-0.08 


lb 


0.02 


0, 


.20 


zb 


0.03 


0, 


.28 


zb 


0.02 


49 


1, 


.73 


± 


0.02 


2, 


.69 


lb 


0.03 


4.43 


lb 


0.02 


9.8 


zb 


0.4 


-0.07 


zb 


0.01 


0, 


.18 


zb 


0.03 


0, 


.25 


zb 


0.02 


50 


1, 


.56 


± 


0.01 


2, 


.70 


lb 


0.02 


4.56 


lb 


0.01 


13.0 


lb 


0.4 


-0.12 


zb 


0.01 


0. 


.15 


zb 


0.01 


0. 


.27 


zb 


0.01 


51 


1, 


.59 


± 


0.02 


2, 


.81 


lb 


0.03 


4.77 


lb 


0.02 


12.0 


lb 


0.4 


-0.04 


zb 


0.01 


0. 


.23 


zb 


0.02 


0. 


.27 


dz 


0.02 


52 


1, 


.61 


± 


0.04 


2, 


.57 


lb 


0.04 


4.48 


lb 


0.03 


12.7 


lb 


0.8 


-0.19 


zb 


0.02 


0. 


.12 


zb 


0.03 


0. 


.31 


zb 


0.02 


53 


1, 


.58 


± 


0.03 


2, 


.65 


lb 


0.03 


4.62 


lb 


0.02 


12.8 


lb 


0.6 


-0.14 


zb 


0.01 


0, 


.16 


zb 


0.03 


0, 


.30 


zb 


0.02 


51 


1. 


,51 


± 


0,01 


2. 


,69 


± 


0,01 


4.80 


± 


0,01 


13,4 


± 


1,0 


-0.12 


± 


0,02 


0, 


.20 


± 


0.0:^ 


0, 


.32 


± 


0.03 



